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ABSTRACT 


Context. Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules 
since they allow important physical properties of the gas and its evolutionary history to be deduced. 

Aims. We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of 
rate coefficients and to respond to the large increase in the number of molecules detected in space since our last release in 2013. 
Methods. We present the sixth release of the UMIST Database for Astrochemistry (UDfA), a major expansion of the gas-phase 
chemistry that describes the synthesis of interstellar and circumstellar molecules. Since our last release, we have undertaken a major 
review of the literature which has increased the number of reactions by over 40% to a total of 8767 and increased the number of species 
by over 55% to 737. We have made a particular attempt to include many of the new species detected in space over the past decade, 
including those from the QUIJOTE and GOTHAM surveys, as well as providing references to the original data sources. 

Results. We use the database to investigate the gas-phase chemistries appropriate to both O-rich and C-rich conditions in TMC-1 and 
to the circumstellar envelope of the C-rich AGB star IRC+10216 and identify successes and failures of gas-phase only models. 
Conclusions. This update is a significant improvement to the UDfA database. For both the dark cloud and C-rich circumstellar 
envelope models, calculations match around 60% of the abundances of observed species to within an order of magnitude. There are 
a number of detected species, however, that are not included in the model either because their gas-phase chemistry is unknown or 
because they are likely formed via surface reactions on icy grains. Future laboratory and theoretical work is needed to include such 


species in reaction networks. 
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1. Introduction 


The first release of the UMIST Database for Astrochemistry 
(UDfA) was made public in 1991 (Millar et al. 1991). It was, 
in main part, motivated by the various astrophysical applications 
— dark clouds, hot cores, circumstellar clouds, novae, and super- 
novae — that were being studied by the UMIST (University of 
Manchester Institute of Science and Technology) Astrochemistry 
Group and the recognition that members of the group should 
have access to one set of reactions and rate coefficients that 
encompassed their various needs. Once we had done this, it 
became clear that releasing it to the wider community would 
encourage the spread of astrochemical modelling as a tool 
both to interpret and to understand and make predictions for 
observations of molecules in space. 

It has been known for some time that gas-phase synthe- 
sis dominates the formation of many important interstellar 
molecules, such as CO, N», HCO*, N5;H*, H5D'*, and the 
unsaturated hydrocarbon chains. For others, however, a detailed 
description of interstellar chemistry has to be augmented through 
reactions within and on the icy mantles of interstellar dust 
grains. Gas-phase reactions do remain, however, the foundation 
on which such ice chemistry rests. It provides the feedstock for 
the grain mantle and can chemically process material that is 
removed from ices. The detections of around 100 new molecules 
in space in the last few years have led to the addition of many 


new reactions and species to gas-phase networks and reflects the 
importance of chemical networks in the interpretation of molec- 
ular line observations. This is particularly true in the current 
paradigm in which *bottom-up' synthesis determines the abun- 
dances of larger species. Readers wishing to keep abreast of 
the latest detections of molecules in space should consult the 
Cologne Database for Molecular Spectroscopy! (Endres et al. 
2016) or the Astrochemyst website? managed by David Woon. 

It is in this context that we present the sixth release of the 
UDfA, RATE22 (previous releases: RATE91 — Millar et al. 1991; 
RATE95 — Millar et al. 1997; RATE99 — Le Teuff et al. 2000; 
RATEO6 — Woodall et al. 2007; RATE12 — McElroy et al. 2013). 
The year suffix refers to the date at which we stopped col- 
lecting new data. Thus, RATE22 contains data published or in 
press up to the end of 2022. Our fifth release (McElroy et al. 
2013) contained 6173 reactions among 467 species and 13 ele- 
ments and was developed in the context of the astronomical 
identification of around 150 molecular species at that time? and 
has proven a very popular source of astrochemical data with 
over 600 references to date. Applications of the data occur to 
sources as diverse as protoplanetary disks (Walsh et al. 2015), 
infrared dark clouds (Entekhabi et al. 2022), brown dwarf disks 


! https: //cdms.astro.uni-koeln.de/ 
? https://www.astrochymist.org/ 
3 http://udfa.net/ 
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Table 1. Code, reaction type, and the number of each reaction type in 
RATEI2 and RATE22. 


Code Reaction type Count 
RI2 R22 
AD Associative detachment 132 142 
CD Collisional dissociation 14 14 
CE Charge exchange 579 . 663 
CP Cosmic-ray proton (CRP) 11 12 
CR Cosmic-ray photon (CRPHOT) 249 428 
DR Dissociative recombination 531 831 
IN Ion-neutral 2589 3447 
MN Mutual neutralisation 981 1501 
NN Neutral-neutral 619 1018 
PH Photoprocess 336 509 
RA Radiative association 92 152 
REA  Radiative electron attachment 24 30 
RR Radiative recombination 16 20 
Total reactions 6173 8767 
Total species 467 737 
DOIs 6990 
URLs 474 
No references 1303 
Measured 1889 
Calculated 1508 
Literature 2996 
Estimated 2374 


Notes. Further information on the method by which the rate coefficient 
has been determined and the source of the data is given in Sect. 2.1. 


(Greenwood et al. 2017), the dynamical and chemical evolu- 
tion of prestellar cores (Priestley et al. 2023), photodissociation 
regions (Róllig & Ossenkopf-Okada 2022), external galaxies 
(Shimonishi et al. 2020), and cometary comae (Heritier et al. 
2017). Its data have also been used in non-astrophysical applica- 
tions, most commonly in plasma physics. 


2. The RATE22 database 
2.1. Description of the data 


Our basic gas-phase ratefile, RATE22^, now contains some 8767 
individual rate coefficients. These correspond to 737 species 
involving 17 elements, increases of over 40% and 55% in reac- 
tions and species, respectively, from RATE12. The additional 
elements are Al, Ar, Ca, and Ti. The basic format is that each 
line of data consists of 18 colon-separated entries: the first two 
are the reaction number and reaction type, defined in Table 1, fol- 
lowed by two reactants and up to four products. The ninth entry 
denotes the number of temperature ranges, NTR, over which the 
rate coefficient is defined, while entries 10—12 give the values 
of o, B and y used to calculate the rate coefficients. Entries 
13-14 give the temperature range over which the rate coeffi- 
cient is defined, entry 15 provides the method by which the rate 
coefficient has been determined (M: measured; C: calculated; 
E: estimated; L: literature). Literature values have generally been 
harvested from reaction networks published within papers or in 


^ https://umistdatabase.net 
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their supplementary materials. It also refers to rate coefficients 
for which we have no information on the method by which they 
have been determined. In some cases, it describes data taken 
from ratefiles sent to us on request by their authors. An exam- 
ple of the latter is the OSU high-temperature network (Harada 
et al. 2010). We note, however, that the difference between esti- 
mated and literature values has been rather eroded over time. Of 
those reactions labelled ‘E’, some 1501 are mutual neutralisa- 
tion reactions, all of which are given the same rate coefficient 
(Harada & Herbst 2008; Loomis et al. 2016). Entry 16 is an esti- 
mate of accuracy, entry 17 gives, where available, some 7464 
Digital Object Identifiers (DOI) or web page (URL) links that 
will take the user to the original data source. As part of this 
update, we have calculated the rate coefficients of 788 ion-dipole 
reactions using the approach described in Sect. 3. Of these, 601 
currently contain no reference label; we shall update these DOI's 
upon publication of this paper which is the source of these rate 
coefficients. Entry 18 lists additional notes in which 'Millar i-d' 
identifies these 788 ion-dipole reactions. A further 78 reactions 
identified by ‘Millar’ refer to rate coefficients calculated or esti- 
mated by us. Of our total set of reactions, only 325 lack both 
a reference and an associated note. Further details on the file 
format are given in McElroy et al. (2013). Table 1 lists the dif- 
ferent reaction types and compares the numbers of each in the 
RATE12 and RATE22 databases. In addition, it gives data on the 
origin and methodology by which the rate coefficients have been 
determined. 


2.2. Calculation of rate coefficients 


A specific rate coefficient is calculated from a, B, and y accord- 
ing to its reaction type. For binary reactions, this corresponds to 
the de Kooij-Arrhenius (KA), or modified Arrhenius, formula: 


TY TY 3-1 
k= a(z) exp( =) ems , (1) 


where T(K) is the gas temperature. For cases in which a single 
rate coefficient is fit by different formulae over several con- 
tiguous temperature ranges (NTR » 1), the entry is extended 
by including the appropriate data as additional colon-separated 
entries that have the same format as entries 10-18. For exam- 
ple, the rate coefficient of the ion-neutral reaction between H ^ 
and CH5, with NTR = 3, has been fitted to three KA formulae, 
one each for the temperature ranges 30-100 K, 100—300 K and 
300—1000 K (Yurtsever et al. 2020). In some complex cases, the 
rate coefficient of a particular reaction is best fitted by a sum 
of KA formulae over the same temperature range. In such cases 
we list the reaction as a series of independent rate coefficients 
each with its specific parameters. An example is the reaction 
between CN and HC3N to form NC4N. A combined experimen- 
tal and theoretical study of this reaction over the temperature 
range 5-400 K by Cheikh Sid Ely et al. (2013) shows that the 
rate coefficient is better fit as the sum of two KA formulae each 
of which is included as a separate line in the ratefile. In total, 
10 reactions have complex rate coefficients described as the sum 
of two or more Arrhenius formulae so that the number of fully 
independent reactions in RATE22 is 8757. 

For one-body reactions, that is those involving photons or 
cosmic ray particles. For cosmic-ray ionisation (CP), the formula 
for evaluation of the rate coefficients becomes 


key s“, Q) 
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where œ is the cosmic-ray ionisation rate, while for those 
involving UV continuum photons (PH) the formula is: 


k = aexp(-yAv) sl, (3) 


where a represents the rate coefficient in the unshielded inter- 
stellar ultraviolet radiation field, Ay is the dust extinction at 
visible wavelengths, and y is the parameter used to take into 
account the increased dust extinction at ultraviolet wavelengths. 
We have not explicitly included the effects of self-shielding in 
the database. Such a process can occur in situations in which 
dissociating photons are absorbed through line rather than con- 
tinuum processes and acts in addition to the extinction caused 
by dust grains. It can be important for molecules including Ho, 
CO, N2, OH, and H2O and depends on the column density 
of the molecule with these species showing almost complete 
shielding once a column density of around 10? cm"? is reached 
(Heays et al. 2017). Numerical approaches often involve the 
use of look-up tables. Details on the self-shielding of H5, CO, 
and N3 are discussed by Sternberg et al. (2014), Visser et al. 
(2009), Li et al. (2013) and Heays et al. (2014). The inclusion of 
self-shielding is particularly important in the atomic to molecu- 
lar transition regions in molecular clouds, in photon-dominated 
regions (PDRs) circumstellar envelopes, and in certain regions 
of protoplanetary disks. 

For cosmic ray-induced photoreactions (CR) the rate coeffi- 
cient becomes: 


T V 
c7 (as) 


where œ is the cosmic-ray ionisation rate, y is the efficiency 
factor as defined in Eq. (8) of Gredel et al. (1989), and w is 
the dust-grain albedo in the far ultraviolet, typically 0.4—0.6 at 
150 nm for particles large compared to the wavelength and close 
to zero for very small grains (Jones et al. 2013). We choose w = 
0.5 in our model calculations. The particular value of y depends 
on the Lyman band photons emitted by H» following collisional 
excitation by the energetic electrons released in cosmic-ray ion- 
isation. Thus the intensity and wavelength dependence of the 
cosmic ray-induced UV flux is very different from interstellar 
photons. Note that our y values are relative to the H» density and 
are different by a factor of two from those listed by Heays et al. 
(2017) to account for the factor 1/(1 — w) in Eq. (4). The cosmic- 
ray ionisation rates listed here are normalised to a total rate 
for electron production from cosmic ray ionisation (primarily 
from H» and He in dark clouds) of Z = 1.36 x 107 s^! (Prasad 
& Huntress 1980). Rates for both direct cosmic ray ionisation 
and cosmic ray-induced photoreactions can be scaled to other 
choices of the ionisation rate, Z, by multiplying the appropriate 
rate coefficients by Z/£o. 

While the absolute rates for CR reactions are different from 
those of PH reactions for reasons mentioned above, we note that 
the number of CR reactions is significantly less than those of 
PH. We have looked in detail at the reasons for this. Three fac- 
tors emerge: (1) ionisation by cosmic-ray photons is negligible 
for some species that readily ionised by the interstellar radiation 
field (see Heays et al. 2017 for examples); (ii) CR destruction of 
molecular ions, particularly where these react rapidly with Ho, 
is ignored; (iii) our literature searches showed that many hydro- 
carbon chain species have several product channels when subject 
to interstellar photodissociation but only one when cosmic-ray- 
induced photons are involved. If we consider photorates for the 
species C,H,, n = 6-11, m = 0-2, there are 61 PH channels 


y -1 
l-w ib 


(4) 
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compared to 21 CR channels (Bettens & Herbst 1995; Wakelam 
et al. 2010; Harada et al. 2010). Since it is not clear how the 
branching ratios might differ between PH and CR reactions, we 
have decided to retain, rather than augment, the original data. 


2.3. Caveats and limitations 


Very recently, Tinacci et al. (2023) used quantum chemistry 
to investigate the thermochemistry of some 5768 reactions and 
the electronic energies of over 500 species contained in the 
KIDA database? (Wakelam et al. 20152) identifying 306 that are 
endothermic. We have searched for each of these in our database 
finding that 53 overlap in terms of reactants and products with 
our list and should be removed. The fact that our number is 
so much smaller than in the KIDA database is due to a num- 
ber of reasons: over 170 of the endothermic reactions identified 
in KIDA, particularly those involving large hydrocarbon chains, 
are not in RATE22; RATE22 does not differentiate between cer- 
tain isomers, for example between 1-C35H and c-C3H nor between 
the cumulene H2C, and diacetylene HC4H, whereas KIDA does; 
energy barriers are already present in the RATE22 reactions; and 
the fact that RATE22 uses lower energy isomers where KIDA 
does not, for example SiOH*, HOSO*, HSCO*, and HCSH* in 
RATE22 versus HSiO*, HSO,*, HOCS*, and H3CS* in KIDA. 

We note too that many of the reactions excluded by Tinacci 
et al. (2023) are ion-neutral reactions with rate coefficients 
measured in low-pressure experiments to be a significant 
fraction of the collisional rate coefficients. In these cases ionic 
products are normally identified by mass spectrometry so that 
a mis-identification of an ionic product may simply be due 
to the fact that the structures identified in the KIDA database 
are not the lowest energy isomer. This is likely to be the case 
where independent measurements give similar products and rate 
coefficients. In other cases, particularly where the branching 
ratio of a channel is small at room temperature, an energy barrier 
may indeed be present. 

As discussed in Sect. 4, the removal of these reactions 
does not make a significant difference to the comparison 
between observed and calculated abundances in both TMC-1 
and IRC+10216. Both the full and the reduced ratefiles are made 
available to the community. 

In addition to the differences between KIDA and RATE22 
discussed above, and noting that there are significant overlaps 
in reactions and rate coefficients, there remain other important 
differences. The KIDA database includes several networks that 
incorporate both gas-phase and grain-surface reactions as well 
as several that describe chemistry in planetary atmospheres. Our 
new database includes reactions involving many of the complex 
organic molecules detected in recent years. It can also be applied 
to the circumstellar chemistries around both C-rich and O-rich 
AGB stars and will, no doubt, continue to be used in the range 
of applications listed in Sect. 1. 

As mentioned in Sect. 1, many situations in astrochemistry 
need to involve surface or bulk chemistry on ice-mantled dust 
grains. The most important surface reaction is the association of 
two H atoms to form H, on a grain surface. This reaction is not 
included explicitly in the RATE22 database. Our software does, 
however, include a term for this process in the ODEs describing 
the formation and destruction of H atoms. That is, it contains a 
term, HLOSS, which accounts for the loss of H atoms through 
the surface formation of H5, the abundance of which is calcu- 
lated, like that of electrons, by a conservation equation. The 
formation of H» on grains has been the subject of many studies 


5 https://kida.astrochem-tools.org/ 
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including the pioneering paper by Hollenbach & Salpeter (1971). 
For example, Cazaux & Tielens (2004) described the process 
by which both physisorption and chemisorption contributed to 
H; formation on bare silicate and carbonaceous surfaces (see 
the correction in Cazaux & Tielens 2010), while Le Bourlot 
et al. (2012) discussed production in PDRs through both the 
Langmuir-Hinshelwood and Eley-Rideal mechanisms and Foley 
et al. (2018) investigated Hz formation on PAH particles. When 
using the database in other codes (that is, not the codes provided 
by us), care must to be taken to include grain-surface formation 
of H, (either implicitly or explicitly). 

Other than H, formation, grain surface chemistry is 
neglected here. It is an essential part of many astrochemical 
applications, particularly in the description of hot cores, hot 
corinos and protoplanetary disks. Many current models do incor- 
porate such reactions but their networks have been developed 
mostly in isolation, that is, specific networks have been designed 
by individuals or groups for a particular purpose; there is no 
central source of 'agreed' data. In particular, a quantitative 
description of surface chemistry needs to include parameters 
such as binding energies, diffusion barriers, reaction products 
and branching ratios, thermal and non-thermal desorption mech- 
anisms, and the nature and composition of the ice, amongst 
others. The method of solving the grain chemistry ODEs is also a 
matter of choice, ranging from the computationally fast, but not 
always applicable, rate equation approach, to slower but more 
accurate Monte Carlo methods. Cuppen et al. (2017) discuss 
these issues in an extensive review article. Experimental and 
theoretical approaches that minimise uncertainties are underway 
but the lack of systematic results has prevented implementa- 
tion of more general rules that might allow surface chemistry 
to be better described. We note, however, that the binding ener- 
gies of some 133 molecules has been published by Ligterink & 
Minissale (2023). Nonetheless, it has long been known experi- 
mentally that binding energies are dependent on the specific site 
at which a molecule sits within the ice (Collings et al. 2004). 
More recent calculations by Bovolenta et al. (2022) and Tinacci 
et al. (2022) have provided the binding energies of over 20 
molecules on model water ice clusters and show that all have a 
broad distribution of binding energies rather than a unique value. 

As discussed above, values of the a, 6, and y parameters 
used to calculate rate coefficients are defined over specific tem- 
perature ranges. In many cases, these are fits to experimental 
data and it is appropriate to give a note of caution to their use 
at temperatures outside these ranges. Róllig (2011) noted that 
the extrapolation to low temperatures of rate coefficients in the 
RATEO6 release could led to divergent and unphysical behaviour 
due to the incorporation of large negative y values. We corrected 
this and other issues noted by Róllig (2011) in our RATEI2 
release and have been careful to ensure that our current rate 
coefficients are not divergent. The choice of whether or not to 
extrapolate a given rate coefficient outside of its given tempera- 
ture range, particularly when determined experimentally, is not 
straightforward and it is unlikely that any one prescription works 
for all reactions. It is known, for example, that the quantum 
tunnelling of H atoms in reactions involving OH can increase 
rate coefficients substantially at low temperatures (Heard 2018). 
Choices also have to be made where rate coefficients are defined 
over two non-contiguous temperature ranges. For this situation 
we recommend that the rate coefficient is interpolated over the 
intermediate range, a solution that has the advantage that it does 
not allow a discontinuity in the value of the rate coefficient. 
Róllig (2011) provides an excellent summary of the issues 
involved in such cases. 
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Itis also possible that certain rate coefficients with large pos- 
itive values of 6 can diverge and be unphysical when extrapolated 
to temperatures much above the upper limit, Ta, of the range over 
which the fit is made (see Shaw et al. 2023 for a discussion of 
this issue in RATE12). This occurs only for NN reactions and 
we have identified 10 such reactions. There is no unique way 
to choose a formulation of these rate coefficients above T, and 
indeed, given the large values of 8 involved in some, it can be 
dangerous to extrapolate to higher temperatures. We therefore 
adopt a simple approach. For these reactions, we set NTR - 2, 
k(T > Ty) = k(T,) and noted this fact in the individual entries. 

Finally, it is worth recalling that the number of reactions that 
describe the chemistry of individual species varies enormously 
from several hundred in some cases, for example CO, H20, C2H, 
and C4H amongst others, to a handful, particularly for many of 
the complex organic molecules and the large, newly detected, 
hydrocarbons. In the former case, one may expect that calculated 
abundances will be less sensitive to the values of rate coeffi- 
cients. Particular circumstances can make this expectation fail. 
An example is diffuse cloud chemistry in which photodissoci- 
ation dominates the loss of neutral species. In this case, their 
abundances should be inversely proportional to their unshielded 
photodissociation rates. In the latter case, the formation routes 
of large molecules are often unknown and have allocated rate 
coefficients that may be highly uncertain. The abundances of 
these species are likely to be almost directly proportional to the 
adopted rate coefficients although complex re-formation path- 
ways may mitigate this behaviour to some extent. We plan a full 
sensitivity analysis of the RATE22 network for both dark clouds 
and circumstellar envelopes in order to identify critical reactions 
for further study. 


3. New and updated reactions and species 


Here we present brief summaries of the chemistry associated 
with new species included since the RATEI2 release as well as 
important updates to reaction rate coefficients. We have made 
a comprehensive literature search for species detected in space 
since our last release, identifying around 100 new molecules. 
We have searched for their formation routes in the literature and 
added destruction through a set of standard loss mechanisms: 
proton transfer reactions with species such as H}, H3O*, HCO*, 
and N5H*, and reactions with He* and C* which tend either to 
break molecular bonds or to build complexity. Charge exchange 
(CE) reactions with H*, and to a lesser extent C*, can also play 
an important role in chemistry. To help augment the formation 
routes of complex species, we have also added some mutual neu- 
tralisation (MN) reactions between common anions observed in 
IRC+10216 and TMC-1 - CN , CN , CN, CH , C6H7, 
and CgH^ - and molecular cations (see Millar et al. 2007, 
Cordiner & Millar 2009 and Walsh et al. 2009). For these, we 
have adopted the rate coefficients suggested by Harada & Herbst 
(2008). Calculated abundances are particularly sensitive when a 
species has only one or a few known formation routes. Since this 
applies to many of the new molecules, we have allowed the disso- 
ciative recombination (DR) of their protonated forms to produce 
smaller fragments as well as the parent neutral. This approach is 
needed in dark clouds to prevent unphysical cycling of a species, 
that is, it prevents the effective loss rate of a molecule from 
being zero. In photon-dominated regions, photodissociation pre- 
vents such recycling. We have included neutral-neutral (NN) 
reactions where we have evidence for these and destruction by 
interstellar and cosmic ray-induced photons. We have searched 
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Table 2. Relabelled species in RATE22. 


RATE12 RATE22 RATE12 RATE22 
H3CO* H5COH* CH;CCH* H5CCCH* 
CH,CO* H5CCO* H3CS* H,CSH* 
C3H,0t HC30H* H3C30* HCCCHOH* 
COOCH,* HCOOCH3* HOCS* OCSH* 
HSiS* SiSH* H5C505* H5COOCH;* 
HSO;* HOSO* C;Hs* CH3C,¢H2* 
H35C;N* HC;NH* H2CoNt HC yNH* 
HC4H C4H2 


for permanent electric dipole moments and provide a list of over 
270 values on our website. We then used these in an approxima- 
tion of the Su-Chesnovich formula (Su & Chesnavich 1982) to 
calculate ion—neutral (IN) rate coefficients: 


k = 3.87 x 10 ?upuc!7(T/300) "2? cm?! , (5) 


where up is the electric dipole moment of the neutral molecule 
in Debye and yp is the reduced mass of the reactants in atomic 
mass units. This approximation is likely correct to within 20% 
for species with dipole moments greater than 1 Debye and is 
well within the uncertainties of astrochemical model calcula- 
tions. Where a specific rate coefficient has been measured at 
300 K, we have used the measured rather than the calculated 
value to scale to other temperatures. 

We have also taken the opportunity to relabel a small number 
of species in RATE12. These are mostly protonated ions and have 
been relabelled to make clear their structural form. Table 2 lists 
these. 

In the following, all fractional abundances are given rel- 
ative to H5 unless otherwise stated. Within each sub-section, 
the discussion of the chemistry is ordered by the mass of the 
molecule. 


3.1. HCCNC, HCCNCH" , HC4NC, HC5NH*, C7N and HC, N 


HCCNC and HCCNCH* were first detected in TMC-1 by 
Kawaguchi et al. (1992) and Agündez et al. (2022a), respec- 
tively, at fractional abundances of 3 x 107!? (Cernicharo et al. 
20202) and 3 x 107? (Agündez et al. 2022a) and an abundance 
ratio HCCNC/HNC; = 5.8 (Cernicharo et al. 2020a). We note 
that Vastel et al. (2018) have detected both HCCNC and HNC; 
in L1544 with an abundance ratio around 10. 

The formation of HCCNC is thought to be due to the DR 
of protonated HC3N. We include five DR channels with rate 
coefficients taken from Wakelam et al. (2010) and adopt the 
same branching ratios for HCCNC and HNC3. Unfortunately, 
the experimental studies of the DR products of DC3ND* do 
not distinguish between the DC3N isomers (Vigren et al. 2012). 
We include proton transfer reactions to form HCCNCH" with 
rate coefficients calculated using Eq. (5) as well as destruction 
involving reactions with C* and He*. 

Isocyanodiacetylene, HC4NC, has been detected in TMC-1 
by Xue et al. (2020) and Cernicharo et al. (2020a) with the latter 
authors determining a fractional abundance of 3 x 107 !!, much 
less than the value of 1.8 x 107? found for HC3N. Protonated 
cyanodiacetylene, HC5NH', has been tentatively detected in 
TMC-1 with a fractional abundance of 7.5 x 10^!! by Marcelino 
et al. (2020). The basic chemistry is taken from Cernicharo 


et al. and, following these authors, we assume that the disso- 
ciative recombination of HC4NH* produces HC4NC with 1% 
efficiency. We have also included formation of HC4NC through 
MN between HC4NC* and the anions listed above. These partic- 
ular processes may be important in both interstellar clouds and 
C-rich AGB outflows such as that around IRC+10216. 

Cernicharo et al. (2023) have detected several transitions of 
the anion CzN in both TMC-1 and IRC+10216 allowing for 
accurate abundance estimates in both sources. It is a new species 
in RATE22 and we have included its formation through the reac- 
tion of H^ with HC7N with a rate coefficient calculated by 
Gianturco et al. (2017) as well as through the radiative electron 
attachment between C7N and e^ (Cordiner & Millar 2009). Loss 
occurs through photodetachment of the electron, reactions with 
H and O atoms as well as in several MN reactions with cations. 

HCN is the largest cyanopolyyne detected to date and was 
identified in TMC-1 by Loomis et al. (2021) with an abundance 
of 10x 10719. It was not included in RATEI2. Loomis et al. 
(2016) provide a detailed chemical model for the synthesis of 
HC,,N as well as several other large carbon-chain molecules, 
including HCzN and HCN, CH3C;N, and CH3CoN, and radi- 
cals such as Cij9N and C; N. We have included close to 100 of 
their reactions in the database to enhance the formation of these 
large species. 


3.2. NC4NH* 


Following the detection in TMC-1 of NCCNH", which is a tracer 
of the unobservable — at least through its rotational line emission 
— cyanogen molecule NCCN, with an abundance 8.6 x 107? by 
Agündez et al. (2015a), a related ion, protonated dicyanoacety- 
lene, NC4NH'*, has been detected in the same source with a 
slightly smaller abundance 1.1 x 107? by Agúndez et al. (20232). 
The chemistry of cyanogen and its related ions has been well 
studied in the laboratory and is described in detail in RATEI2. 
Its larger counterpart NC4N is likely to form through exothermic 
NN reactions involving C3N and both HCN and HNC, studied 
in detail by Petrie & Osamura (2004). A more efficient route 
though is the reaction of CN with HC3N which has been stud- 
ied experimentally over the temperature range 22-296 K by 
Cheikh Sid Ely et al. (2013). These authors also fit their data 
to higher temperature measurements and fit a double Arrhenius 
form from 5-400 K to the rate coefficient which has a value in 
excess of 107!? cm? s^! at 10 K. 


3.3. Propargyl and related species: CH2CCH and CH2C3N 


Propargyl, CH2CCH, has been detected in TMC-1 with a frac- 
tional abundance of 107? by Agúndez et al. (2021a, 2022b) 
making it one of the most abundant hydrocarbon radicals yet 
detected. It was included in RATE12 but its ion-neutral reactions 
were erroneously given a 7-7! dependence. Since its dipole 
moment is small, 0.15 D, this dependence has been removed. 

3-cyanopropargyl, CH2C3N, has been detected in TMC-1 by 
Cabezas et al. (2021b) with a fractional abundance of 1.6 x 107!!. 
This molecule is newly included in the RATE22 database and we 
follow the approach discussed by Cabezas et al. with formation 
driven by four NN reactions: 


C + CH;CHCN > CH>C3N + H (6) 
C; +CH3CN > CH;C;N +H (7) 
CN + CH;CCH 2 CH>C3N + H (8) 
CH; + CjN > CH>C3N + H, (9) 
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with rate coefficients as suggested by Cabezas et al. Loss is 
assumed to be via the standard set of ion-neutral reactions with 
the DR of protonated CH,C3N assumed to be divided equally 
between the product channels CH2C3N + H and CH3 + HC3N. 


3.4. Methyl cyanopolyynes and their isomers 
3.4.1. CH3C3N, HCCCH2CN and H2CCCHCN 


Two isomers of methyl cyanoacetylene, CH3C3N, propargyl 
cyanide, HCCCH2CN, and cyanoallene, H: CCCHCN, have been 
detected in TMC-1 by Marcelino et al. (2021). All three species 
have very large dipole moments and we have calculated ion- 
neutral rate coefficients for new chemistry involving these 
species using Eq. (5). We have, in the absence of evidence, 
assumed that the DR of the protonated forms of both propar- 
gyl cyanide and cyanoallene results in their parent neutral and 
propargyl, CH2CCH, with a 50:50 branching ratio. Formation 
of all three species occurs primarily through neutral-neutral 
reactions involving CN: with propyne, CH3CCH, and allene, 
CH5CCH», to form CH3C3N (Abeysekera et al. 2015; Carty et al. 
2001); with allene to form propargyl cyanide at a branching ratio 
of 0.1 (Balucani et al. 2002), and with both propyne and allene 
to form cyanoallene, with branching ratios again taken from 
Balucani et al. (2002). 


3.4.2. CH3C5N and H2CCCHC3N 


CH3Cs5N was detected in TMC-1 by Snyder et al. (2006). 
Fuentetaja et al. (2022b) used their most recent data to calcu- 
late its fractional abundance as 9.5 x 107 ?, about 8 times smaller 
than the value given by Snyder et al. (2006), arguing that their 
analysis fits better to the high-J transitions that they observe. 
Its isomer, cyanoacetyleneallene H;CCCHC3N, was detected 
in TMC-1 via the line stacking method by Shingledecker 
et al. (2021) who derive a preferred value of its abundance as 
2 x 107!!. Fuentetaja et al. (2022b) subsequently confirmed this 
identification through the detection of 20 rotational transitions 
and derived an abundance of 1.2 x 107!!. 

RATE12 assumed that the radiative association between 
CH;* and HC3N leads to CH3C5NH* followed by DR to 
form CH3C5N (Herbst & Leung 1989). For its formation, we 
add the detailed pathways and rate coefficients suggested by 
Shingledecker et al. (2021) and based on the laboratory study 
by Fournier (2014) who studied the reaction of C3N with allene 
and propyne down to 24 K. The allene reactions are assumed 
to produce H2CCCHC3N with 50% efficiency. Shingledecker 
et al. assume that there are five product channels, with equal 
branching ratios, of the propyne reaction, including the for- 
mation of both CH43C3N and H;CCCHC3N. Two of these, 
however, are species not contained in the database, namely 
CH3CC and HCCCH2C3N. We thus substitute these species 
as their equivalents, propargyl CH2CCH, and H,CCCHC3N. 
In addition to these two reactions involving C3N, we include 
another neutral—neutral formation reaction, that between CN and 
acetylene allene, H;CCCHCCH. Since the dipole moment of 
H5CCCHOCSN is large, 5.15 D (Cabezas, priv. comm.), we have 
also included the usual ion-neutral reactions with enhanced rate 
coefficients. 


3.4.3. CH3C7N and CH3CgN 


CH3C7N was discovered in TMC-1 by Siebert et al. (2022a) 
who determined an abundance of 8.6 x 107. Its chemistry in 
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RATE12 is parallel to that of CH3C5N but with the radiative 
association now involving HC7N rather than HC5N (Herbst & 
Leung 1989). 

CH3CoN is newly included in the database and is as yet unde- 
tected in TMC-1. Its formation is assumed to be through the 
radiative association of CH3* and HCoN with a rate coefficient 
assumed to be equal to that of the RA forming CH3C;NH'. It 
has a large dipole moment, 6.5 D, and we included the usual set 
of ion-neutral reactions with rate coefficients determined from 
Eq. (5). CH3C9NH* is formed in proton transfer reactions as 
well as in the reaction of N with Cj9Hs* and is lost via DR 
and by reaction with C atoms to form H5C;; NH* (Loomis et al. 
2016). For both species we have added some new ion-neutral 
loss reactions and calculated their rate coefficients using Eq. (5). 


3.5. CH2CHCCH 


Vinyl acetylene was included in the RATEI2 release but its 
chemistry was extremely limited with destruction only via 
interstellar and cosmic-ray-induced photons. Here, we have 
added formation through the measured C2H + CH, reaction 
(Bouwman et al. 2012), destruction through IN and NN reac- 
tions involving C,H to form ortho-benzyne, c-CgH4 (Zhang 
et al. 2011), and CN to form, as suggested by Lee et al. 
(2021b), the more complex and newly-detected interstellar 
species vinylcyanoacetylene, H;CCHC3N, and cyanovinylacety- 
lene, HCCCHCHCN, the chemistry of which are discussed 
further in Sect. 3.18. 


3.6. Methyl polyacetylenes and their isomers 
3.6.1. CH3C4H and H2CCCHCCH 


Allenyl acetylene, H;?CCCHCCH, was detected in TMC-1 
by Cernicharo et al. (2021f) with a fractional abundance of 
1.2x 1072, essentially identical to that of its isomer methyl 
diacetylene, CH3C4H. Here we adopt the formation reactions 
suggested by Cernicharo et al., namely reactions of C,H with 
CH5,CCH, and CH3CCH, and of C2 with CH3CHCH» which 
have been measured at low temperatures (Carty et al. 2001; 
Daugey et al. 2008). Product branching ratios for these two reac- 
tions have been determined by Goulay et al. (2007). The dipole 
moment of allenyl acetylene is small, 0.63 D, so its ion—neutral 
reactions do not have a T~'/? dependence. 


3.6.2. CH3CgH and HeCCCHC4H 


Allenyl diacetylene, H;CCCHCAH, an isomer of CH3C&H, has 
been detected in TMC-1 with an abundance 2.2 x 107!" by 
Fuentetaja et al. (2022b) who also determine the abundance of 
CH3C¢H to be 7x107!!. We adopt the same neutral-neutral 
formation routes as Fuentetaja et al. which are based on the low- 
temperature experiments of Berteloite et al. (2010). The RATE12 
database contains many measured ion-neutral reactions that pro- 
duce the C7H;* ion. Unfortunately, information on the isomeric 
forms of this ion produced in the experiments is unknown. Here 
we assume that the isomeric forms are split equally between 
CH;C6H2* and H5CCCHC4H»*. Note that in RATEI2 proto- 
nated CH3C&H was labelled as C7Hs*. We explicitly differentiate 
between the protonated forms of the two isomers in RATE22. 


3.6.3. CH3CgH 


Siebert et al. (2022a) searched unsuccessfully for CH3CgH in 
TMC-1 deriving an upper limit to its abundance of 9.8 x 107!°. 
We have included it in RATE22 assuming that it is one of the 
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products of the dissociative recombination of Co9Hs* formed in 
the radiative association of CsH3* and C4H» (Herbst & Leung 
1989). The DR rate coefficient and its branching ratios are set 
equal to those of C7H5*. The C9Hs* ion is also lost through reac- 
tion with N atoms, leading eventually to HCoN (Loomis et al. 
2016). 


3.7. Molecules containing the C—-S bond 
3.71. HCS and HSC 


The first detections of the HCS radical and its metastable isomer 
HSC were made by Agúndez et al. (2018b) in the dark cloud 
L483. Cernicharo et al. (2021g), in their survey of S-bearing 
molecules in TMC-1, give an abundance of 5.5 x 107!° for HCS 
and an upper limit of 1.3 x 107!! for HSC. 

We have added formation of HCS via the fast reactions 
C + H2S (Deeyamulla & Husain 2006) and S with propargyl 
as well as some DR channels in the recombination of larger 
S-bearing ions. We note that RATEI2 included photoionisation 
but not photodissociation of HCS. We now include dissociation 
from both interstellar and cosmic ray-induced photons, setting 
the rates for these equal to those for HCO. The gas-phase syn- 
thesis of HSC is essentially unknown (see Agündez et al. 2018b) 
and we do not include it here. 


3.7.2. HCCS, HC5S*, H2CCS, HC3S, H2CCCS, C48, HC4S 
and C5S 


A number of new S-bearing hydrocarbons were discovered in 
TMC-1 by Cernicharo et al. (2021g), who discuss their chemistry 
in detail. We have adopted their approach to the chemistry, taking 
many of the reactions and their rate coefficients from Vidal et al. 
(2017) in their pioneering study of complex S-chemistry. Since 
all species have dipole moments greater than 1 D, we have re- 
calculated all ion-neutral reactions using Eq. (5), adding proton 
transfer and IN reactions involving abundant ions where these 
are missing. 

HC5S* has been detected in TMC-1 by Cabezas et al. 
(2022b) at an abundance of 1.1 x 10719. It appears in the RATE12 
database but we have added a small number of reactions here. As 
C5S has a very large proton affinity, 869.6 kJ mol", larger than 
that of NH3 and smaller than that of C3S, we have added some 
additional proton transfer reactions for the formation of HC2S* 
and included proton transfer to C3S as a loss mechanism. 

The HC,S molecule has been detected in TMC-1 with an 
abundance of 9.5 x 107? by Fuentetaja et al. (2022a). We include 
the primary formation routes suggested by them, notably the 
reactions of S and S* with C4H5* and C4H3, respectively, and 
C with H2CCCS. To these we have added formation through S 
and C4H». Loss reactions include those with N atoms and charge 
exchange and proton transfer reactions with rate coefficients 
calculated using Eq. (5). 

The thio-carbon chains, C4S and C3S, which have been 
detected in TMC-1 with fractional abundances of 3.8 x 107? 
for C4S and 5.0x 107? for CsS by Cernicharo et al. (2021), 
have large dipole moments and therefore the IN reactions involv- 
ing both species have a 7-!? dependence at low temperature. 
C4S was included in RATE12 with its formation dominated by 
MN reactions between anions and C4S*. We have recalculated 
all IN rate coefficients using Eq. (5) and added some neutral- 
neutral formation routes including the reactions between S and 
C4H and between C and HC5S (Vidal et al. 2017) as well as 
C and H5CCCS. C58 is a new species in the database and we 


have followed closely the approach for its formation and destruc- 
tion suggested by Cernicharo et al. and Vidal et al., although we 
have again calculated IN rate coefficients using the ion-dipole 
approach. 


3.7.3. NCS, HNCS and HSCN 


Isothiocyanic acid, HNCS, and its metastable isomer thiocyanic 
acid, HSCN, were detected first in the Sgr B2 molecular cloud by 
Frerking et al. (1979) and Halfen et al. (2009), respectively. Sub- 
sequently, Adande et al. (2010) detected both species in TMC-1 
while Cernicharo et al. (2021g) detected NCS in the same 
object and gave abundances for all three species, 7.8 x 10^ !!, 
3.8 x 107!! and 5.8 x 10^!! , respectively. The chemistry of these 
three molecules are heavily interlinked and has been discussed in 
some detail by Gronowski & Kotos (2014), who focused primar- 
ily on the routes forming HNCSH* and H2NCS*, and Vidal et al. 
(2017). We adopt their pathways but for consistency recalculate 
IN rate coefficients using Eq. (5). 


3.7.4. HCSC?H and HCSCN 


Propynethial, HCSC?H, and cyanothioformaldehyde, HCSCN, 
have been detected in TMC-1 by Cernicharo et al. (20211) 
with abundances of 3.2 x 107!! and 1.3 x 10716, respectively. For 
HCSC2H we have included the formation reactions suggested 
by Cernicharo et al., namely S with CH?CCH and C2H with 
H5CS, with rate coefficients of 107! cm? s~!. We have included 
ion-neutral destruction with He* and C* as well as proton trans- 
fer reactions with the main proton carriers in dark clouds, H3*, 
HCO*, N;H*, and H;0*, followed by dissociative recombination 
with two equal product channels to HCSCoH + H and C2H3 + 
HCS. To enhance its formation, we have also included a number 
of MN reactions as discussed in Sect. 1. 

The formation of HCSCN is treated in a similar manner to 
that of HCSC?H, that is, we adopt two formation reactions S + 
CH5CN and CN + H5CS. These reactions, and those mentioned 
above for HCSC2H, are analogues of the O and H2CO reactions 
that produce the interstellar molecules HCCCHO and HCOCN 
(see Sect. 3.9). 


3.7.5. CH3SH and CoHsSH 


Methyl mercaptan or methanethiol, CH3SH, was first detected in 
the interstellar medium by Linke et al. (1979) and has been seen 
in a wide variety of massive star-forming regions as well as in 
the low-mass star-forming system IRAS 16293-2422 (Majumdar 
et al. 2016). It has not been detected in the classical cold dark 
clouds TMC-1 and L483. Two recent papers (el Akel et al. 2022; 
Bouscasse et al. 2022) have performed extensive line surveys 
of two complex sources containing both hot and cold gas to 
unravel sulphur chemistry. Bouscasse at el. find CH3SH only 
in hot gas surrounding G328.2551-0.5321 while el Akal et al. 
find it in the cold envelope of Cyg X-N12 as well as in the 
hot gas. The formation of CH3SH in the gas-phase appears to 
be difficult with no efficient routes identified that involve either 
ion-neutral or neutral-neutral reactions, leading to the conclu- 
sion that its chemistry is dominated by grain surface reactions 
(see, for example, Majumdar et al. 2016). 


One measured ion-neutral reaction that forms the 
C-S bond is: 
S+ + CH3CHO — CH3SH* + CO, (10) 


although the branching ratio to this channel is only 0.2 (Decker 
et al. 2000). We assume that the CH3SH* ion is converted to 
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CH3SH through a series of MN reactions with anions. Loss 
of CH3SH occurs in ion-neutral reactions. Where there is no 
experimental evidence, we have adopted the branching ratios 
suggested by Majumdar et al. although we have been more con- 
servative in including reactions and products that are analogous 
to reactions involving CH3OH. 

Ethyl mercaptan, C2Hs5SH, has been detected in its gauche 
(lowest energy) form toward the galactic centre by Rodríguez- 
Almeida et al. (2021a) who discuss its formation via possible 
grain-surface reactions (see also Lamberts 2018). Due to the lack 
of feasible gas-phase formation routes, we have not included this 
species in the database. 


3.8. The hydrocarbon chains — CsH*, C5H», CgHo, 
HCCCHCCC, CioH and CoH- 


CsH* was detected in TMC-1 by Cernicharo et al. (20222) at an 
abundance of 8.8 x 107?. CsH* has many formation routes but 
relatively few destruction reactions in RATE12. We have added 
proton transfer reactions involving Cs with No2H* and H3O* 
to the former and included MN loss reactions with anions to 
the latter. 

The linear cumulene pentatetraenylidene, CsH2, has been 
detected with one of the lowest abundances seen in TMC-1, of 
1.35 x 107? (Cabezas et al. 2021d). Its chemistry was covered 
extensively in RATE12 but, as its dipole moment is 5.9 D, we 
have updated all IN rate coefficients involving this species using 
Eq. (5). 

Cabezas et al. (2021d) also reported the abundances of the 
ortho- and para-cumulene forms of C6H2 with an ortho-para 
ratio of 3 and a total abundance equal to 8x 107? in TMC-1 
while Fuentetaja et al. (2022a) identified several emission lines 
in the same source from its isomer HCCCHCCC with an abun- 
dance of 1.3 x 107!!. We do not include the latter species in the 
network. 

The largest pure hydrocarbon chains detected to date are 
CjoH and Cj9H both of which have been observed, the former 
a tentative detection, by Remijan et al. (2023) in TMC-1 with 
abundances of 2.0 x 10^!! and 4.0 x 10^ !!, respectively. As noted 
by Remijan et al. (2023), this is the first neutral-anion pair in 
which the anion has a larger abundance than the neutral, a result 
difficult to explain with our current, and still limited, under- 
standing of anion chemistry. These two species were already 
included in RATE12 but their chemistry was incomplete since, 
unlike smaller hydrocarbon chains, we did not allow their for- 
mation from larger chains. We have rectified this by including 
CuH and Cu H^ in our species set and including the reaction 
O + CH” to form CjgH- + CO. The formation of a smaller 
hydrocarbon anion has been found in the measurements of O 
with CH , CH , and CH. by Eichelberger et al. (2007) who 
also note that the reactions become more exothermic as the chain 
length increases. Although we fit the observed abundance of 
CjoH well, we find, in common with Remijan et al. (2023), that 
we underproduce Cj )H™ (see Sect. 4.1.1). 


3.9. HCCCHO, c-C3H20, H2CCCO, HCOCN and 
CH2CHCHO 


Propynal, HCCCHO, was detected in TMC-1 by Irvine et al. 
(1988) and subsequently observed by Cernicharo et al. (20211) 
who derived a fractional abundance of 1.5 x 1079. The chem- 
istry of its isomers, c-C3H2O (cyclopropenone), with an abun- 
dance of 5.4x107!? (Loison et al. 2016), and H2CCCO 
(propadienone), which has an upper limit of 1.1 x 107"! 
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(Cernicharo et al. 2020b) have been discussed in some detail by 
Loison et al. Here, we include only the chemistry of HCCCHO 
and, in particular, the two reactions considered by Cernicharo 
et al. (20211): 


O + CH?CCH — HCCCHO + H 
C,H + H5CO — HCCCHO + H, 


(41) 
(12) 
with rate coefficients of 107!? cm? s~!. Ion-neutral rate coeffi- 
cients are calculated using Eq. (5) since the dipole moment of 
propynal is 2.78 D. 

Remijan et al. (2008) made the first detection of formyl 
cyanide (or cyanoformaldehyde), HCOCN, in Sgr B2(N) and 
suggested that it was formed through the reaction of CN with 
H2CO, a suggestion strengthened by the quantum theoretical 
study of Tonolo et al. (2020). More recently, however, West et al. 
(2023) studied the reaction experimentally over the range 32- 
103 K and performed high-level quantum calculations to find that 
the dominant product channel is HCN + HCO. They found that 
the pathway to HCOCN formation is hindered by a submerged 
energy barrier. Based on a combination of experimental and the- 
oretical approaches, they produced a fit to the rate coefficient 
over the temperature range 6—1500 K. In RATE22, the synthesis 
of HCOCN is dominated by two NN reactions: 


O + CH?CN  HCOCN + H 
CN + CH3CHO  HCOCN + CH3, 


(13) 
(14) 


with rate coefficients equal to 10-!° (Cernicharo et al. 2021i) and 
6.5 x 107? cm? s7! (Tonolo et al. 2020). In TMC-1, Cernicharo 
et al. find the fractional abundance of HCOCN is 3.5 x 107!!. 

The lowest energy conformer of propenal, CH; CHCHO, was 
originally discovered by Hollis et al. (2004) in the giant molec- 
ular cloud Sgr B2(N) and subsequently detected in TMC-1 by 
Agúndez et al. (2021b) with an abundance of 2.2 x 107!!. It is 
formed through NN reactions (Tsang 1991; Goulay et al. 2012; 
Xie et al. 2005): 


O + C3Hs > CH;CHCHO + H (15) 
CH + CH;CHO > CH;CHCHO + H (16) 
H5CO + C,H; > CH;CHCHO + H. (17) 


We also include its formation through MN between the com- 
mon anions and CH?CHCHO*, formed in charge exchange with 
H* and C*. 


3.10. HəNC and H2CN 


H2CN was included in RATE12 but its chemistry was limited and 
dominated by neutral-neutral reactions taken from Smith et al. 
(2004). Since its dipole moment is 2.54 D, we have added our 
usual suite of ion-neutral reactions with rate coefficients calcu- 
lated by Eq. (5). H2NC, the highest energy isomer of H2CN, was 
discovered in the dark clouds L483 and B1-b by Cabezas et al. 
(20212) although they did not detect it nor H2CN in TMC-1, 
with upper limits to their fractional abundances of 3.2 x 10-!! 
and 4.8 x 107!!, respectively. Although the energy difference is 
29.9 kJ mol!, the abundance ratio of these two molecules is 
close to 1 in both cold clouds. We have adopted a similar chem- 
istry for H,NC, adding the measured reaction between C and 
NH; (Bourgalais et al. 2015) to form both isomers with equal 
efficiency. Neutral atom destruction of H2NC is assumed to 
proceed at the same rate as that of HCN. 
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3.11. C,O and HCO, n = 4-9 


In recent years, there has been a number of successful (and 
unsuccessful) observational searches in TMC-1 for long carbon- 
chain molecules terminated by an oxygen atom (McGuire et al. 
2017a; Cordiner et al. 2017; Cernicharo et al. 2021b). Between 
them, these papers give fractional abundances for species CO 
to C7O and HCO to HC;O. 

Two major mechanisms have been suggested to form these 
molecules. The first is based on experimental work on the ternary 
association reactions between CO and unsaturated hydrocarbon 
ions from which are deduced two-body radiative association 
rate coefficients at 10 K to form HC,,,;0*, HC,,,OH* ,and 
H2Cn+1OH* ions, n = 2-6 (Adams et al. 1989). We have adjusted 
their pre-exponential factor and adopted a (T/300) ?? depen- 
dence so that the rates can be used at other temperatures. Such a 
dependence is based on the modified thermal approach (Harada 
et al. 2010) where the exponent depends on the sum of the rota- 
tional degrees of freedom of the reactants. In this case we assume 
that the hydrocarbon ions are non-linear. If instead, the dom- 
inant collisional ion is linear, then the dependence should be 
(T/300)-??. This would mean that we overestimate the rate con- 
stant at 10 K by a factor of 5.5. In addition, we have extended the 
set of reactions to include the C9O species. Since the 10 K rates 
are already close to their collisional values for n = 5, we assume 
that the radiative association rates for these larger species are 
equal to that for the CSH* + CO association. 

Once formed these ions can then undergo DR to form the 
range of neutrals observed. It should be clear, however, that the 
actual abundances of the (H)C,O species are particularly sen- 
sitive to the branching ratios for these DR processes. In the 
absence of any evidence, we have chosen a procedure that is 
easy to implement (and to change in the light of future informa- 
tion) — each ion formed through radiative association dissociates 
into two channels with equal branching ratios. Using C7O as an 
example: 


HC70* + e7 > Ce + CO 0.5 (18) 
HC;O* +e —> CO +H 0.5 (19) 
HC7OH* + e > C70 + H5 0.5 (20) 
HC;OH* + e. > HC;O + H 0.5 (21) 
H;C;OH* +e” 2 C70 +H +H 0.5 (22) 
H5C;OH* + e > HC;O + H5 0.5. (23) 


The second route to these long chain molecules has been 
proposed by Cordiner & Charnley (2012) and is based on the 
laboratory measurements of Eichelberger et al. (2007). These 
form the chain species through reactions involving O atoms and 
hydrocarbon anions, such as: 


O + C6H7 — HC6O + e7 
O +CH —> CO +C;H7, 


(24) 
(25) 


again with equal branching ratios. 

In both of these scenarios, a certain simplicity of imple- 
mentation is adopted. The real situation is likely to be more 
complex with other possible product channels and branching 
ratios so our approach, in the absence of further investigation, 
likely overestimates the gas-phase formation of these neutral 
chain molecules. 

The C40 species all have large dipole moments (Moazzen- 
Ahmadi & Zebretto 1995); those of HC,O are smaller 
(Mohamed et al. 2005) but still larger than 1 D. We use Eq. (5) 
to generate the rate coefficients for ion-neutral reactions. 


3.12. CH3CHCHCN 


The trans- and cis- forms of crotononitrile, CH3CHCHCN, 
are two of five cyano derivatives of propene discovered in 
TMC-1 by Cernicharo et al. (2022c), the others being methacry- 
lonitrile, CH;CH(CH3)CN, and gauche- and cis-allyl cyanide, 
CH2CHCH2CN. All species have very similar abundances, 
—10-!, indicating that they are likely to form in the same 
manner. Cernicharo et al. note that the reaction of CN with 
propene, CH;CHCH2, has been measured to be fast down to 
23 K (Morales et al. 2010) with products that could include 
the above species as well as vinyl cyanide, CH2CHCN, with a 
branching ratio that greatly favours vinyl cyanide. 

Since the detailed chemistry of these cyano derivatives and 
the interplay between them is not understood, here we only 
include crotononitrile as a proxy for the others and allocate it 
a branching ratio of 30%, with vinyl cyanide allocated 70% in 
the CN-CH3CHCH) reaction. In addition to photodissociation, 
we include the usual ion-neutral loss reactions. As noted by 
Cernicharo et al., this formation route is now a very inefficient 
mechanism for forming these cyano derivatives. This is a result 
of a very much decreased production rate of propene follow- 
ing the study by Lin et al. (2013) who showed that the radiative 
associations of both C3;H3* and C3H;* with H» possess energy 
barriers and do not proceed at low temperatures. The result is 
a precipitous decrease in the abundance of C3H7* and hence of 
propene calculated under dark cloud conditions, well below that 
observed, and as a result, leads to very low abundances for these 
cyano derivatives (see Sect. 4.1.1). 


3.13. Cyanates, their isomers and related species 
3.13.1. HNCO and its isomers 


Isocyanic acid, HNCO, has three metastable isomers whose 
chemistry has been discussed in detail by Quan et al. (2010). 
Here we updated ion-neutral rate coefficients using Eq. (5) and 
added proton transfer reactions with HCO*, N2H*, and H30*. 


3.13.2. OCN, HoNCO* and H20CN* 


These molecules have not been observed as yet in TMC-1 but the 
first two were detected in L483 by Marcelino et al. (2018). OCN 
has a fairly extensive chemistry in the RATE12 database but the 
only proton transfer reaction included was with H3* which erro- 
neously had a 7T-!? dependence. We have corrected this and 
included a wider range of such reactions. We have also added 
a new formation reaction CN + O, for which Sims et al. (1992) 
measured the rate coefficient from 26 to 295 K. At higher tem- 
peratures we adopt the rate coefficient recommended by Baulch 
et al. (1994). 

The chemistry of H2NCO* and H,OCN% is very limited in 
RATE12 where they are formed only by proton transfer of H3* 
with HNCO and HOCN respectively, and by the reactions of H2 
with HCNO* and HOCN', respectively. Marcelino et al. have 
shown that the H» reactions are, however, endothermic. We have 
therefore removed the H» reactions and added additional proton 
transfer reactions. Furthermore, the only loss routes for the two 
ions in RATEI2 are via DR; we have added MN reactions as 
a result. 


3.13.3. CHgNCO, HOCH2CN and C2H5NCO 


The first detection of interstellar methyl isocyanate, 
CH3NCO, was made in Sgr B2(N) by Halfen et al. (2015). 
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Majumdar et al. (2018) have discussed its gas and grain 
chemistry. Their list of reactions contains only one gas-phase 
formation pathway, the reaction of HNCO and CH3. Their cal- 
culation of the rate coefficient, however, has an energy barrier 
of 8040 K making it unimportant in most interstellar sources, 
with the possible exception of shocked gas. Their modelled 
abundance of CH3NCO therefore depends completely on grain 
surface reactions. Subsequently, Gorai et al. (2020) revised the 
value of the rate coefficient to a temperature-independent value 
of L0 x 107? cm? s*!, whereas Quénard et al. (2018) adopt a 
value of 5.0 x 107!! cm? s^!. Since this is the single gas-phase 
production reaction for CH3;NCO, its calculated abundance is 
directly proportional to the value of the rate coefficient. Here we 
adopt the value of Quénard et al. to maximise its abundance but 
note that this will be overestimated by a factor of 50 if Gorai et 
al.’s value is correct. We have taken the bulk of our additional 
gas-phase reactions from Majumdar et al., updating them using 
Eq. (5) where appropriate. 

Glycolonitrile, HOCH2CN, an isomer of CH3NCO, has been 
detected in two interstellar sources IRAS 16293-2422B (Zeng 
et al. 2019) and Serpens SMM 1 (Ligterink et al. 2021). There are, 
however, no gas-phase routes identified to form this molecule 
and we do not include it in the ratefile. 

Finally, we note in passing that ethyl isocyanate, C2H5NCO, 
has been detected in the galactic centre cloud G+0.693-0.027 by 
Rodriguez-Almeida et al. (2021b). There is no gas-phase chem- 
istry available for the formation of this molecule and it is not 
included in the database. 


3.14. NH2CHO 


Formamide, NH2CHO, is observed in many interstellar clouds 
although not in TMC-1, where its upper limit is 5x 107? 
(Cernicharo et al. 2020b), and was not included in RATE12, con- 
sistent with the subsequent findings by Redondo et al. (2014b,a) 
that potential gas-phase ion-neutral syntheses had energy bar- 
riers despite being exothermic. The molecule was added to the 
current database following the calculation by Barone et al. (2015) 
that the reaction between NH; and H5CO was efficient in making 
NH2CHO. However, a very recent experimental and theoreti- 
cal study of this reaction down to 34 K shows that the channel 
producing NH;CHO has an energy barrier of over 1800 K and 
therefore cannot act as a production route in interstellar clouds 
(Douglas et al. 2022). There appears to be no route to NH;?CHO* 
thereby ruling out MN as a low-temperature route to formamide. 
We do, however, leave this species in the database since the high 
temperatures that need to be reached to overcome the energy bar- 
rier can be realised in hot, inner circumstellar envelopes (CSEs) 
and in post-shock gas. 


3.15. HCOOCHs, CH3COOH, HOCH2CHO and (CHOH)2 


Acetic acid, CH3COOH, glycolaldehyde HOCH2CHO, and 1,2- 
ethenediol, (CHOH), hereafter AA, GA, and ED, respectively, 
are isomers of the common interstellar molecule methyl formate, 
HCOOCHs, hereafter MF, the only one of the four included 
in RATE12. Its formation there, though, is very inefficient, 
essentially due to one slow ion-neutral reaction, that between 
CH3OH»* and H5CO to produce protonated MF. The RATE12 
model produces a maximum abundance of less than 107! under 
TMC-1 conditions (Sect. 4.1.1). MF is the only one of these 
four species to be detected in TMC-1 with an abundance of 
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1.1 x 10770 (Agündez et al. 2021b). The others are observed in 
molecular clouds associated with star forming regions. 

Ascenzi et al. (2019) have measured the rate coefficient of 
He* with MF and re-evaluated other proton transfer reactions. 
We have followed their suggestions but used Eq. (5) to cal- 
culate total rate coefficients except where measured values are 
available. 

Gas-phase formation of AA and GA, as currently under- 
stood, is an inefficient process and most models that attempt to 
reproduce their abundances in interstellar clouds tend to rely on 
grain surface chemistry. Skouteris et al. (2018) discuss this in 
some detail. Their gas-phase scheme begins with the hydrogen 
abstraction reactions of Cl and OH with ethanol: 


Cl + CH5OH — CH3CHOH + HCl 
Cl + C.Hs0H — CH2CH2OH + HCI. 


(26) 
(27) 


These radicals then react with O atoms — CH3;CHOH to form 
AA and CH;CH5,OH to form GA, amongst other products. 
Skouteris et al. calculate rate coefficients and branching ratios 
for all these reactions which we adopt here. As usual, we have 
used Eq. (5) to calculate the rate coefficients of the ion—neutral 
loss reactions. We have taken DR rate coefficients for both 
ionised and protonated AA and GA from Walsh (priv. comm.) 
with some adjustments to branching ratios for channels involv- 
ing products CH3CO and COOH that are not contained in 
this release. 

(Z)-1,2-ethenediol, which we write in the database as 
HOCHCHOH, has recently been detected in the Galactic Centre 
cloud G+0.693-0.027 by Rivilla et al. (2022). They find a low 
rotational excitation temperature of about 8.6 K, and an abun- 
dance ratio with respect to glycolaldehyde of 0.19. We include 
the two gas-phase routes suggested by Rivilla et al.: 


H2CO + CHOH 5 HOCHCHOH + H (28) 
CH;CHOH + OH — HOCHCHOH + H, (29) 
with — temperature-independent rate coefficients 2x 


10719 cm? s-!. 


3.16. CH3NH, CH2NH» and CH3NH» 


Methylamine, CH3NHz2, was one of the earliest detected inter- 
stellar molecules having been observed in the Giant Molecular 
Clouds Sgr B2 and Orion A by Kaifu et al. (1974) and Fourikis 
et al. (1974), respectively. It was not included in earlier versions 
of the database. Its gas-phase formation is highly uncertain with 
the dominant pathway initiated by the radiative association of 
CH;* and NH; to form CH3NH;3* (Herbst 1985) followed by 
dissociative recombination. It has recently been shown experi- 
mentally that it has a fast loss reaction with OH at 22 K in which 
the radicals CH2NH» and CH3NH are produced (González et al. 
2022) with a rate coefficient some 20 times larger than that at 
300 K. Puzzarini et al. (2020) have made a theoretical study of 
the H-atom abstraction reaction of CH3NH», with CN and we 
adopt their rate coefficients and branching ratios to CH; NH» and 
CH3NH. These two radicals have also been added to RATE22. 
We note that Schmid et al. (2022) have recently obtained a 
laboratory rotational spectrum of protonated methylamine and 
searched for it in several cores in Sgr B2, determining an upper 
limit of ~107!°. Finally, we note that Wlodek & Bohme (1988) 
have studied the reaction of Si* with CH3NH» showing that it 
forms SiNH* and hence SiN following DR. 
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3.17. Cyclic hydrocarbons — c-Cs5H, c-C3HCCH, c-C5He, 
C-CgHa4, c-C6Hs, c-C5H4CCHp, C-C5H5CCH, C-C5H5CN, 
CgH5CCH, c-CgHsCN, c-CgHg and C-CgH;CN 


c-CsH, written as C4CCH to distinguish it from linear C5H in 
the database, has been detected with an abundance of 9 x 10? 
in TMC-1 (Cabezas et al. 2022a). They suggest two formation 
reactions, one the reaction of C2H with c-C3H, which retains the 
triangular C3 structure in c-C5H. The second is the reaction of 
C atoms with diacetylene, a reaction they suggest may produce 
both linear and cyclic C5H, although they note that production 
of the latter is close to thermo-neutral and may not proceed at 
low temperatures. We have therefore adopted their suggestion 
and taken a branching ratio of 10% for c-Cs5H. In addition to 
destruction by N, C and O atoms, we have added our usual suite 
of dipole-enhanced IN reactions. 

Ethynyl cyclopropenylidenene, c-C4HCCH, written simply 
as C3HCCH in the database, and cyclopentadiene, c-C5Hg, and 
written CsHg, were detected in TMC-1 with fractional abun- 
dances of 3.1 x 107!! and 1.2 x 107°, respectively, by Cernicharo 
et al. (2021e,c). Since the basic structure of the former species 
is a C3-triangle, its chemistry is thought to be dominated by 
the reaction of C;H with cyclopropenylidene, c-C3H», given that 
both reactants are abundant in TMC-1 and that other reactions 
between C2H and hydrocarbon radicals are known to be fast at 
low temperature. c-C3HCCH has a large dipole moment, 3.54 D, 
so again we have used Eq. (5) to calculate ion-neutral rate coeffi- 
cients. Theoretical and/or experimental studies on reactions that 
might lead to this molecule are clearly needed. 

The gas-phase chemistry of cyclopentadiene, c-C5Hg, is 
also likely to be incomplete. Its neutral-neutral formation is 
entirely dependent on one reaction between CH and butadi- 
ene, CH;CHCHCH, (Cernicharo et al. 2021c). The dissociative 
recombination of CsH7* can also produce it but, other than pro- 
ton transfer reactions of c-C5sH forming CsH;*, which are cyclic 
in nature and not therefore true formation routes to C5H6, this 
ion is only produced in slow ion-molecule reactions involving 
C3H4* and C3Hs*. In general, the unreactivity of large carbon- 
bearing species with H5 makes the synthesis of such molecules 
with more than 3—4 hydrogen atoms, and in particular C5H;*, 
inefficient. 

Cernicharo et al. (2021d) detected ortho-benzyne, hereafter 
C6H4, with a fractional abundance of 5.0x 1071! in TMC-1. 
They provide an exhaustive discussion of a total of 14 possi- 
ble formation routes, many of which possess energy barriers. 
Here, we include the four major neutral-neutral routes identi- 
fied by Cernicharo et al. The exothermic, barrierless reaction 
between C2H and vinyl acetylene, CH;CHCCH, has been mea- 
sured by Zhang et al. (2011) with an estimated rate coefficient of 
8 x 1071! cm? s~!. The other three syntheses involve the reactions 
of C;H4, C3H, and C3H, with C4H, CH3CCH, and CH,CCH, 
respectively. We have also included the reaction between C3H 
and CH5CCH, with all rate coefficients set to 107!?cm?s7!. 
Loss reactions include charge exchange with H* and Ct and 
proton transfer reactions. The phenyl radical, hereafter C6H5, 
has not yet been detected in the interstellar medium but plays 
an important role in the chemistry of the benzene-like ring 
molecules. 

Fulvenallene, CSH4CCH», was detected with an abundance 
of 2.7 x 1071? in TMC-1 by Cernicharo et al. (2022b). They sug- 
gested that it was formed as one of the products in the reaction 
of CoH with cyclopentadiene, c-C5He, the other two channels 
producing 1- and 2-ethynyl cyclopentadiene, c-C4H5CCH. Since 
both propene and c-C5Hg play an important role in the synthesis 
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of fulvenallene we have added a number of ion-molecule reac- 
tions, taken from the compilations by Anicich & Huntress (1986) 
and Anicich (1993), aimed at increasing the abundances of these 
species. These include: 


CoH? + CH3CH3 = C4H; + CH, (30) 
C3Hz + CoH4 > C5H; + H5 (31) 
CH;CCH'* + C2H4 > CsHz + H2. (32) 


Cernicharo et al. (2022b) suggest that the radiative asso- 
ciation reaction between I-C3H3*, denoted CH»CCH* in the 
database, and C2H4 forms C5H;* with a very large rate coeffi- 
cient. Although this product ion was detected in the laboratory, 
the measurement was performed under the relatively high pres- 
sure conditions within an ion cyclotron resonance experiment, 
so that the rate coefficient determined is more likely to be the 
two-body equivalent of a three-body process rather than that of 
radiative association. 

The two isomers l- and 2-ethynyl cyclopentadiene, 
CsHsCCH, have been detected in TMC-1 by Cernicharo et al. 
(2021e) with abundances of 1.4 x 107!° and 2.0 x 107!°, respec- 
tively. We do not differentiate between these two forms in the 
network. In addition to the reaction between C2H and c-C5Hg 
mentioned above, we have also included formation via the dis- 
sociative recombination of the C7H7* ion which forms primarily 
from the reaction of CH;CCH* with benzene (Smyth et al. 1982) 
as well as in ion-neutral reactions between c-CgHs* and various 
hydrocarbons (Ausloos et al. 1989). 

The detection of 1-cyanocyclopentadiene, CsHsCN, in 
TMC-1 by McCarthy et al. (2021) was followed quickly by the 
identification of 2-cyanocyclopentadiene (Lee et al. 2021a). Lee 
et al. analysed their velocity-stacked spectra and used MCMC 
modelling to derive total fractional abundances of 8.3 x 107! 
and 1.9 x 10^!!, respectively. We ignore the difference in struc- 
ture between these two isomers in the database. Lee et al. suggest 
that the molecule is formed by the reaction of CN with cyclopen- 
tadiene, c-C5Hg, in analogy with the formation of benzonitrile 
from benzene, while Cernicharo et al. (2021e) suggest forma- 
tion via the reaction between CN and CH; CHCHCH;. Here we 
take the rate coefficient for the former reaction to be the same as 
that for benzonitrile formation. Given their large dipole moments 
(Lee et al. 2021a), proton transfer reactions and other ion-neutral 
rate coefficients are calculated using Eq. (5). 

Cernicharo et al. (2021e) tentatively detected ethynyl ben- 
zene, c-CgHsCCH, in TMC-1 at an abundance 2.5 x 10-9. They 
suggested two formation routes, reactions of C,H with CgHg and 
C&H5CN, where the first of these has been measured to have a 
fast reaction down to 103 K (Goulay & Leone 2006). Although 
not applicable in the case of dark clouds, we have added the mea- 
sured high-temperature formation reaction between C2H2 and 
CeHs to the network and adopted temperature-independent rate 
coefficients for its destruction. 

Benzonitrile, CcH5CN was first detected in TMC-1 by 
McGuire et al. (2018). Subsequently, Burkhardt et al. (2021b) 
detected it in four other cold sources and derived a total abun- 
dance of 1.6 x 107? in TMC-1, using an MCMC analysis to fit 
their data to four components in their line-of-sight. Burkhardt 
et al. present a model for its interstellar synthesis. It is readily 
formed in the reaction of CN with benzene, C6H6, a reaction 
that is measured to be fast down to 15 K (Cooke et al. 2020). 
Both its proton affinity and dipole moment are relatively large 
and we therefore allow proton-transfer reactions with H3*, N2H*, 
HCO*, and H30* to produce CeH5CNH" with rate coefficients 
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calculated using Eq. (5). We note that proton transfer reactions, 
and therefore CcH5CNH', are not included in the Burkhardt et 
al. chemical model. We assume that its DR leads to benzonitrile 
and phenyl, C6H5, with equal branching ratios. 

Indene, c-CoHg, was detected in TMC-1 by both Cernicharo 
et al. (2021c) and Burkhardt et al. (2021a). The latter authors 
have described a pathway to indene based on the work of 
Doddipatla et al. (2021) and involving the reaction of CH with 
styrene, c-CgH5C2H3. Their calculated abundances fall several 
orders of magnitude below that observed. We have not included 
indene in the database. 

Sita et al. (2022) have recently searched for the isomers 
of cyanoindene, C>H7CN in TMC-1, detecting 2-cyanoindene 
at a fractional abundance of 2.1 x 10^!!. The formation of this 
molecule likely occurs in a manner similar to that of benzonitrile, 
that is, through a fast reaction between CN and indene. Such a 
reaction results in an abundance ratio of cyanoindene to indene 
close to that observed but gas-phase synthesis, as currently 
understood, is not able to reproduce the observed abundance by 
several orders of magnitude. Cyanoindene is not included here. 

The model failure to reproduce the large observed abun- 
dances of the benzene-related cyclic hydrocarbons is ultimately 
due to a lack of efficient, barrierless pathways to C6H6, mainly 
due to the fact that the fast pathway to propene formation is no 
longer feasible as discussed in Sect. 3.12 (Lin et al. 2013). The 
identification of such pathways, or alternative routes to benzene, 
is a key issue for astrochemistry research. 


3.18. HoCCHC3N and HCCCHCHCN 


Vinylcyanoacetylene, H;CCHC3N, and its isomer cyanoviny- 
lacetylene, HCCCHCHCN, the latter in its trans-(E) conformer, 
have been detected in TMC-1 by Lee et al. (2021b). They cal- 
culate abundances of 2 x 10^!! for the former and 3 x 107!! for 
the latter as well as an upper limit of 2 x 10^!! for the trans-(Z) 
conformer of cyanovinylacetylene. The neutral-neutral reactions 
that form these species are taken from Lee et al. to which we add 
our usual suite of ion-neutral and photoreactions. 


3.19. CH»CHOH 


Vinyl alcohol, CH CHOH, an isomer of acetaldehyde, 
CH3CHO, has been detected in TMC-1 at an abundance of 
2.5 x 107? by Agúndez et al. (2021b), very close to the abun- 
dance, 3.5 x 107!9, of acetaldehyde (Cernicharo et al. 2020b). In 
general, experimental studies of neutral-neutral reactions cannot 
differentiate between isomeric products. As a result, it is possi- 
ble that those NN reactions that are thought to produce CH3CHO 
could also produce CH?CHOH. The essential equality of the iso- 
meric abundances suggests that both form by the same processes. 
RATE12 already contains several neutral-neutral reactions that 
form CH3CHO so here we make the simple assumption that both 
isomers are produced by these reactions with equal branching 
ratios. We have made the same assumption for the DR of the 
protonated forms of both CH;CHOH and CH3CHO. 


3.20. HCCN, HCCO and HC4N 


Thirty years separate the discovery of HCCN in the outer enve- 
lope of IRC+10216 (Guélin & Cernicharo 1991) and in TMC-1 
(Cernicharo et al. 2021a). Its synthesis has been discussed by 
Osamura & Petrie (2004) and, with particular reference to 
Titan chemistry, by Loison et al. (2015). Its formation is via 
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neutral-neutral reactions: 


CH + HCN 5 HCCN +H (33) 
CH + HNC 5 HCCN +H (34) 
C + HCN 5 HCCN +H, (35) 


with loss through reactions with H, N and CH; and with ions. 

HCCO was first detected in dark clouds by Agündez et al. 
(2015b) and reported in TMC-1 by Cernicharo et al. (2020b, 
2021b) with a fractional abundance of 7.7 x 107!!. Its chem- 
istry was discussed by Wakelam et al. (2015b) who proposed 
both gas-phase and grain-surface syntheses. In the gas phase its 
major formation reactions are those of C;H with OH and O2. For 
the former, we adopt the rate coefficient suggested by Wakelam 
et al. For the latter, we include the four product channels in the 
compilation by Baulch et al. (2005) with that to HCCO given a 
branching ratio of 2096. Loss occurs through H and O atoms 
(Baulch et al. 2005) together with ion-neutral reactions. We 
follow the recommendation by Wakelam et al. (2015b) for the 
products and branching ratios for the dissociative recombination 
of H5CCO*, replacing those contained in RATEI2 which were 
taken from Prasad & Huntress (1980). HCCO was not contained 
in RATEI2 but its protonated form, CH,CO*, was included. 
For consistency, we have relabelled CH2CO* as H5CCO* in 
RATE22. 

HC,N was detected in TMC-1 by Cernicharo et al. (2021a) 
with an abundance of 3.7 x 107!!. We assume that it is formed in 
a similar fashion to HCCN, that is through reactions of CH and 
both HC3N and its isomer HNC;. 


CH +HC3N 2 HC4N + H (36) 
CH + HNC; 2 HC4N + H (37) 
CN +C3H> 2 HCuN+H (38) 
CN + HCCC > HC4N +H (39) 
CN + C2H3 > HC4N +H, (40) 


with rate coefficients taken from Loison et al. (2015). We calcu- 
late loss rates in IN reactions with common ions and also include 
loss with N atoms. 


3.21. CH3CI 


Methyl chloride was detected in the proto-binary source IRAS 
16293-2422B by Fayolle et al. (2017). The chemistry of this 
molecule and related species has been discussed in detail by 
Acharyya & Herbst (2017) and we adopt their scheme for 
RATE22. 


3.22. PH3 


Phosphine was not included in the RATE12 release as no route to 
PH,* could be found. As the PH; reactions with H» are endother- 
mic, formation of the P-H bond proceeds via the RA of P* and 
PH* with H». However, PH3* does not react with H to form 
PH,*. Here we circumvent this issue by incorporating MN reac- 
tions of PH3* with anions to form PH3 directly. Its loss via 
ion-neutral reactions is taken from Charnley & Millar (1994) 
with its photodissociation rate from Sil et al. (2021). 


3.23. ArH* 


ArH* was detected by Barlow et al. (2013) in the Crab Nebula 
and subsequently in absorption toward several Galactic molec- 
ular clouds by Schilke et al. (2014). Argon has an unusual 
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interstellar chemistry in that its proton affinity is lower than that 
of H) while its ionisation potential is larger than that of atomic 
hydrogen meaning that Ar cannot be ionised by the interstellar 
UV field nor can ArH* be formed through the proton transfer 
between Ar and H3*. Instead, ArH* is formed by the reaction 
of Ar* with H5 where the ion is the result of cosmic ray ion- 
isation of Ar and charge exchange between He* and Ar (Babb 
& McLaughlin 2018). We have adopted the chemistry and rate 
coefficients proposed by Schilke et al. augmented by the labo- 
ratory measurements of the rate coefficients of ArH* with H» 
and CO (Villinger et al. 1982). Its DR rate coefficient has been 
calculated by Abdoulanziz et al. (2018) who show that it is van- 
ishingly small for vibrationally cold ArH* and we neglect it here. 
Schilke et al. show that that ArH* is a tracer of gas with a low 
H; fraction. 


3.24. Circumstellar chemistry 


In the following, we describe the chemistries of elements 
that, due to their large depletion in interstellar clouds, are 
predominantly found in molecular form in the circumstellar 
environments. 


3.24.1. MgO and MgOH 


MgO and MgOH are formed through NN reactions in either the 
hot inner winds of O-rich AGB stars or following grain disrup- 
tion in post-shock gas, rate coefficients taken from the theoretical 
calculations by Decin et al. (2018). In general these reactions 
have large energy barriers. MgO also forms in the radiative 
association of Mg and O (Bai et al. 2021a) but the rate coef- 
ficient is too small to produce it in any significant amount. Its 
photodissociation rate is taken from Bai et al. (2021b). 


3.24.2. MgNC, MgCN and HMgNC 


These Mg-bearing species have all been detected in the outer 
envelope of the carbon-rich AGB star, IRC+10216 — MgNC 
(Guélin et al. 1993), MgCN (Ziurys et al. 1995), and HMgNC 
(Cabezas et al. 2013) — indicating that a relatively low- 
temperature chemistry produces these species. Cabezas et al. 
calculate column densities of 1.3 x 10? (MgNC), 74x 10!! 
(MgCN), and 6 x 10!! (HMgNC) cm^ to show their abundance 
ratios scale roughly as MgNC:MgCN:HMgNC = 15:1:0.8. The 
first two are likely formed in the dissociative recombination of 
MgNCH" formed by the radiative association of Mg* and HCN 
(Petrie & Dunbar 2000). We used the Levenberg-Marquardt 
algorithm to fit their tabulated data over the temperature range 
10—100 K by the expression 1.96 x 107! (T/300)7!® cm? s^! and 
added the association of Mg* with HNC which is slightly faster 
than that with HCN. Cabezas et al. also suggest that some 
HMgNC could be formed with a branching ratio of about 196 in 
the dissociative recombination of larger Mg-cyanopolyyne ions 
and we include this possibility (see Sect. 3.24.3). Mutual neu- 
tralisation between the protonated versions of these molecules 
with the observed anions, listed in Sect. 1, are also included. 
For MN reactions involving MgNCH* we have assumed a 50:50 
branching ratio to MgNC and MgCN. 

We have included the standard set of IN reactions and ion- 
dipole rate coefficients for the loss of these species. For the case 
of IRC+10216, however, photodissociation is likely to be a more 
important loss process in the external envelope and one that can 
determine both column densities and the radial extent of these 
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Table 3. Fit parameters for the radiative association of Mg* with 
cyanopolyynes over the range 10-300 K. 


Neutral @ p y 
HCN 196x10-7  —149 0.0 
HNC 371x107  -148 0.0 
HC3N 7715x107^ -1.94 10.79 
HC3N 190x107? -2.24 33.6 
HCN  398x10?  -0.705 0.0 
HCN 6.92x10-? -0.48 0.0 


Notes. The fits for HCN and HNC are over the range 10-100 K. 


species. We have allocated uniform (and highly uncertain) rates 
of 10719e-! 7^v s-! for all these species. 


3.24.3. MgC3N, HMgCsN, MgCsN and MgC;N 


The magnesium cyanopolyynes, MgC3N and MgCsN have both 
been detected in the outer envelope of IRC+10216 with col- 
umn densities of 9.3 x 10? and 4.7 x 10? cm", respectively 
(Cernicharo et al. 2019; Pardo et al. 2021). The formation of 
these species is thought to proceed from the radiative association 
of Mg* with the HC5,,; N cyanopolyynes. Their rate coefficients 
have been tabulated by Dunbar & Petrie (2002) over the range 
10-300 K. Our fits to these are shown in Table 3. For com- 
pleteness, we include the rate coefficient for HC9N in this table 
although MgCoN is not included in the ratefile. On the other 
hand, MgC7N, which is not yet detected, is included here. 

The large ions created via radiative association will disso- 
ciatively recombine with electrons and could, in principle, have 
several product channels. Here, we have followed the scaling pro- 
cedures adopted by Pardo et al. (2021) with a total rate coefficient 
of 3x 1077 (T/300)-!? cm? s*!. As an example, the dissocia- 
tive recombination of MgCSNH* leads to MgNC, MgCsN, 
and HMgNC with values of the rate coefficients at 300 K of 
7.80 x 1075, 2.22 x 1077, and 2.22 x 107? cm? s^, respectively. In 
addition, we have included MN reactions between the observed 
anions and Mg-bearing cations. The Mg-carbon chain molecules 
have very large dipole moments, in excess of 6 D, and there- 
fore very large rate coefficients for ion—neutral reactions at low 
temperature. As for the smaller Mg-CN species, their abundance 
and extent in the outer envelope is likely to be determined by 
photodissociation; we adopt the same photodissociation rates as 
those for the smaller Mg-bearing species. In all cases we assume 
that the product is Mg and the appropriate C5, ,N radical. 

Recently Cabezas et al. (2023) detected HMgC3N in 
IRC+10216 with a column density of 3 x 10? cm". Its forma- 
tion is uncertain but may occur as a product channel in the DR of 
the larger Mg-cyanopolyyne ions. It is not included in RATE22. 


3.24.4. MgC2, MgC2H, MgC4H, MgCeH and MgCgH 


The magnesium acetylide chains MgC;H, MgC4H, and MgC;H 
have all been detected in the outer envelope of IRC+10216 
(Agündez et al. 2014; Cernicharo et al. 2019; Pardo et al. 2021) 
with column densities of 2 x 107, 2.2 x 10!3, and 2 x 10? cm?, 
respectively, and abundance ratios of 1:11:10 in order of increas- 
ing size (Pardo et al. 2021). More recently, MgC; has been 
detected in the same region with a column density of 10? cm~? 
(Changala et al. 2022). Their formation routes are, as for the 
magnesium cyanopolyynes, thought to be dominated by radiative 
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Table 4. Fit parameters for the radiative association of Mg* with 
polyynes over the range 10-300 K. 


Neutral œ p y 
C2H3 9.40x1075 —147 0.0 
C4H2 2.57x107!5 -1.955 12.0 
C6H2 119x107? -2.21 2121 
CH2 341x107 -1.78 33.29 
CoH? 2.28x107? -0.60 10.48 


Notes. The fit for C,H) is over the range 10-100 K. 


association, in this case between Mg* and the polyacetylenes. 
These rate coefficients are tabulated by Petrie & Dunbar (2000) 
and Dunbar & Petrie (2002) and we have again found a best fit to 
them (Table 4). The ions so produced are assumed to recombine 
in a manner in which 80% of the recombinations lose either a 
hydrogen or a magnesium atom, for example, dissociative recom- 
bination of MgC4H»* leads to either Mg + C4H3 or MgC4H + 
H with equal branching ratios, with the remaining 20% divided 
equally between MgC, and MgC;H. Of these species, the most 
difficult to match calculated and observed column densities is 
MgC4H for which the value of the RA rate coefficient is more 
than three orders of magnitude smaller than those for the larger 
hydrocarbons. 

In common with the magnesium cyanopolyynes, photodis- 
sociation plays an important role in setting the abundance and 
radial extent of these molecules. We assume the same rate coef- 
ficients as in Sect. 3.24.2 and assume that the sole product is Mg 
and the appropriate C5, H radical. 


3.24.5. Aluminium chemistry 


Seventeen Al-bearing molecules, including AIO, AIOH, AIF, 
and AICI which are detected in AGB stars (Decin et al. 2017; 
Danilovich et al. 2021), are included in this release. Given the 
fact that along with Ti and Ca, aluminium is one of the most 
depleted elements in interstellar clouds, our primary purpose is 
to include a network that describes as best as possible the syn- 
thesis of observed Al-bearing molecules in hot gas in the inner 
circumstellar envelopes of these stars. As such, it is dominated 
by nearly 100 neutral-neutral reactions and contains molecules 
in size up to 5 atoms (Al,03) and in mass up to 132 amu (AICI;). 
Some of these, particularly the oxides, have measured rate coef- 
ficients but the majority are theoretical estimates which, because 
of the high-temperature environment in which these species are 
found, have both forward and backward determinations of rate 
coefficients. The main sources of information on these reactions 
comes are from Swihart et al. (2003) and Gobrecht et al. (2022). 
We note that the latter authors, who are interested in how Al-O 
clusters can grow to dust grains, describe routes to molecules 
much larger than included in the database. 

There is no gas-phase formation for AIH; included in 
RATE22. We retain it for those occasions in which it may be 
chosen as a parent molecule in the calculation of Al chemistry 
in the inner regions of O-rich AGB stars. If not included as a 
parent, its abundance is always calculated to be zero. 


3.24.6. Calcium chemistry 


Eleven NN reactions are used to describe CaO and CaOH for- 
mation at high temperature, appropriate for conditions in the 
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Table 5. Fit parameters for the radiative association of Na* with 
cyanopolyynes over the range 10—300 K. 


Neutral a« p y 
HCN 2.32x107!8 -1.50 0.0 
HNC 112x107  -150 0.0 
HCjN 349x105  -196 11.25 
HCSN 243x107 -2.65 30.09 
HCN .3148x1079 -2.41 35.11 
HCN  2.65x10?  -124 1718 


Notes. The fits for HCN and HNC are over the range 10-100 K. 


inner envelopes of O-rich AGB stars. Rate coefficients are a 
mix of theoretical, taken mainly from Decin et al. (2018), and 
experimental values. 


3.24.7. Titanium chemistry 


TiO and TiO; are detected species in the circumstellar envelopes 
of O-rich AGB stars (Decin et al. 2020; Wallstróm et al. 2024). 
Laboratory data exists for the synthesis of both through a mix- 
ture of IN and NN chemistry. Rate coefficients for the latter 
are mostly experimental and taken from Campbell & McClean 
(1993) and Plane (2013). Ti is reactive with several common 
interstellar and circumstellar oxides to form TiO which can be 
converted to TiO, through reaction with OH and O5. Some of 
these reactions are known to be efficient at room temperature 
and, since both have large dipole moments, we have included 
proton transfer with rate coefficients calculated using Eq. (5) 
which, following DR, lead to loss mechanisms for both. A low- 
temperature interstellar synthesis of these species is therefore 
possible though this will be limited in effect since Ti is one of 
the most heavily depleted elements in interstellar gas. 


3.24.8. Sodium chemistry 


Petrie & Dunbar (2000) and Dunbar & Petrie (2002) give rate 
coefficients for the RA of Na* with HCN and the cyanopolyynes 
up to HCoN for temperatures up to 100 K and 300 K, respec- 
tively. As for the corresponding data for Mg*, we have fitted their 
tabulated values to the usual modified Arrhenius formula. NaCN 
is a T-shaped molecule and has been observed in IRC+10216, 
originally by Turner et al. (1994) and mapped with ALMA by 
Quintana-Lacaci et al. (2017) who show that, in common with 
CH3CN (Agündez et al. 2015c) and HC3N (Siebert et al. 2022b), 
it has emission close to the star and extends out to a distance of 
a few times 10? cm. 

To estimate the abundance of NaCN, which has a dipole 
moment of 8.85 D, we have included the usual set of ion-dipole 
rate coefficients and adopted an unshielded photodissociation 
rate of 1071? s71. The branching ratios of the DR reactions of 
the NaCzn+ıNH* ions (n = 0-4) are unknown; here we adopt a 
50:50 ratio between the Na and NaCN channels. 

Cabezas et al. (2023) detected NaC3N in the outer enve- 
lope of IRC+10216 with a column density of 1.2 x 10!! cm~?, 
about a factor of eight less than that of MgC3N and a detection 
aided by its extremely large dipole moment, 12.9 D. Its formation 
routes are unknown but, as with the Mg cyanopolyynes, they are 
likely to involve the dissociative recombination of the products 
of the radiative association between Na* and the cyanopolyynes. 
Table 5 shows our fits to the rate coefficients calculated by 
Dunbar & Petrie (2002) for the radiative associations of Na* with 
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Table 6. Detected species, listed by number of atoms, missing from the RATE22 database. 


2/3  SiP KCl SiCN 
KCN CaNC HSC 

4 HOOH CNCN HNCN 

5 CNCHO HCOSH HNCNH 
MgCsN* 

6 CH3NC c-C4H50 CH;CNH 
MgC4H* HMgC3N HOCOOH 

7 CH20CH 2 HOCH2CN HCCCHNH 

8 NH5CH;CN CH3CHNH CH35SiCH; 
MeCoH* 

9 CH3CONH2 C2H5SH CH3NHCHO? 

10 C,H;sCHO (CH;0H); CH3;0CH,0H 
CH2C(CH3)CN C2H5NH2 CoH5NCO 

11 C2H5OCHO CH3COOCH; NH2CH2CH20H 

>12 C34H;CN C;H5OCH; c-C}yH7CN 


SiCSi AINC FeCN 
NCO HSO 

HONO NCCP? 

NH2OH NH2CN NaC3N 
SiH3CN NHCHCN NCCHNH 
H5CCCN 

c-C;H40 MgC;N* 

NH5CONH, HCCCHCCC CH3CHCO 
HOCHCHCHO 

CH,CHCH;0 — CH;CCHC;N CH;CHCH;CN 
C,HsCHO NH5C(O)CH;0H 

C4H5CN CH3COCH20H 

c-CoHg c-CyoH7CN C5H;0H 


Notes. Tentative identifications are labeled by *?'. 


the cyanopolyynes. Below 30 K, the rate coefficients for Na* are, 
however, more than an order of magnitude smaller than those of 
Mg* with HC3N and HCSN indicating that efficient formation 
of NaC3N may depend on either a larger branching ratio or rate 
coefficent in the DR of NaC3NH', or on preferential production 
via DR of the larger RA products. We note that NaC3N is not yet 
included in the database. 

Sodium chloride has been detected most frequently in the 
CSEs of both C-rich and O-rich AGB stars (Cernicharo & 
Guélin 1987; Decin et al. 2016), in the post-AGB object CRL 
2688 (Highberger et al. 2003) and in one early-type region, 
the protostellar disk around Orion SrcI (Ginsburg et al. 2019). 
The formation of NaCl is very poorly understood. Simsová- 
Zámecníková et al. (2022) have calculated the rate coefficient 
of the direct formation through the RA between Na and Cl 
over the temperature range 30-750 K. We have also included 
formation via the reaction of Na and HCI for which we adopt 
a fit over 591—966 K to the two experimental data sets in the 
NIST database. Petrie (1996) has suggested that NaCl can react 
with CN to produce NaCN and we have included this reaction 
with a (very uncertain) rate coefficient of 1071? cm? s7!. NaCl 
has a large dipole moment, 9.0 D, so ion-dipole rate coefficients 
are adopted, with photodissociation and cosmic ray-induced 
photodissociation rates taken from Heays et al. (2017). 

Finally, it is worth pointing out that Na has a very low 
ionization potential. However, Na* is unreactive with many com- 
mon neutral species which makes it an important source of free 
electrons in dark clouds. 


3.25. Other detected molecules 


As noted previously, there are a number of molecules detected 
in recent years that are not included in RATE22. The main 
reasons for such exclusion are: (1) that gas-phase formation 
routes have not yet been identified for them; (ii) their for- 
mation, while known, depends on the presence of another 
species whose origin is unknown; and (iii) grain surface chem- 
istry likely dominates their syntheses. In addition to those 
molecules mentioned previously, these include ethylene oxide, 
also known as oxirane, c-C2H4O, isocyanogen, CNCN, detected 
in TMC-1 by Agúndez et al. (2018a), nitrous acid, HONO, 


3-hydroxypropenal, HOCHCHCHO, also known as malonalde- 
hyde, detected in IRAS 16293-2422B (Coutens et al. 2019, 
2022), and methoxymethanol, CH3;0CH2OH, discovered in 
NGC63341] MM1 (McGuire et al. 2017b). Missing species, 
ignoring conformers and the fullerene molecules, are shown 
in Table 6. 


4. Model results 


In this section, we describe our model calculations for conditions 
appropriate to the dark cloud TMC-1 and to the circumstellar 
envelope of the carbon-rich AGB star IRC+10216. In both cases, 
we discuss the results for two ratefiles — our standard file and the 
reduced file in which the 53 endothermic reactions identified by 
Tinacci et al. (2023) are excluded (Sect. 2.3) and compare these 
to the extensive observational data for both objects. 


4.1. Dark cloud TMC-1 


We choose typical dark cloud parameters; n(H2) = 10*cm™?, a 


kinetic temperature T = 104 K, a visual extinction of 10 mag., 
and a cosmic-ray ionisation rate of 1.3 x 10 s^!. 


4.1.1. O-rich TMC-1 model 


We have calculated a time-dependent model using the ‘low 
metal' initial elemental abundances given in Table 7. These 
O-rich (C/O less than 1) abundances are the same as those in 
McElroy et al. (2013) with the exception of S for which we adopt 
the value of 1.5 x 1076 as suggested by Fuente et al. (2023) in 
their recent modelling of the chemistry of dark clouds in Taurus 
and Perseus. 

Figures 1 and 2 show the time evolution of various N-bearing 
molecules in the O-rich, RATE22 and RATE12 models. There 
are significant differences between the results that are, in major 
part, driven by the inclusion of the reaction of C + NH; in 
RATE22. This reaction was studied experimentally from 50 K 
to 300 K by Hickson et al. (2015) and shown to be fast at low 
temperatures. Since atomic carbon is the dominant form of this 
element before it is processed into CO, the abundance of NH; is 
depressed at times up to about 5 x 10? yr. For the cyanopolyynes, 
the situation is somewhat different. Here, their peak abundances 
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Table 7. Initial abundances relative to total H nuclei, ny. 


Species — nj/ngy) Species i nj/ng 
Hə 0.5 Na 2.0(-09) 
H 5.0(-05) Mg 77.0(—09) 
He 0.00 Si 8.0(-09) 
C 14(-04) P 3.0(-09) 
N 7.5(-05) S 1.5(-06) 
O 3.2(-04) Cl 4.0(—09) 
F 2.0008) Fe 3.0(-09) 


Notes. “'a(b) = a x 10^. 
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Fig. 1. Time evolution of various N-bearing molecules in the dark cloud, 
O-rich model for both RATE22 (solid lines) and RATE12 (dashed lines). 
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Fig. 2. Time evolution of various cyanopolyynes in the dark cloud, 
O-rich model for both RATE22 (solid lines) and RATE12 (dashed lines). 


in RATE22 are reached at slightly earlier times than in RATE12 
but the peak abundances are fairly similar. The major differences 
now occur at long times, greater than 10° yr. 

We have searched for a ‘best time’ fit between calculated, 
fc, and observed, fo, fractional abundances for those 134 species 
for which we have 'retrieved' abundances, that is, we neglect 
those species for which only upper limits are available. To do 
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Fig. 3. Best-fit, D, and reduced best-fit, Dred, average difference 
between observed and calculated species for the full ratefile and 
O-rich initial conditions. The curve labelled ‘NR’ indicates the num- 
ber of species removed in calculating the global average, Dred, and N is 
the number of species whose abundances lie within an order of magni- 
tude of those observed. 


Table 8. Fitting information on the comparison between observed and 
calculated abundances for the O-rich case. 


Time (yr) D Dyed NR N 

6.3x10^  2.03(2.14) 1.38(L44) 16(18) 52(45) 
1.6 x 10° 1.89(1.95)  L41(1.46) 9(10)  55(53) 
5.0x10?  2.51(2.52) L73(172) 21(Q2) 35(35) 


Notes. Values in parentheses are values when 53 reactions identified by 
Tinacci et al. (2023) are set equal to zero. See text for further details. 


so, we calculate the average modulus of the difference, D, in the 
logarithms of these between these two abundances: 


D= 


X([logf. — log fo]?)'”. (41) 


obs 


Since we know that gas-phase chemistry is not appropriate for 
certain species - CH3OH is a well-known example — we also cal- 
culated a ‘reduced’ difference, Drea, in which we remove species 
with D > 4, that is, those whose calculated abundance differs 
by more than 4 orders of magnitude from that observed. Figure 3 
shows values of D, Dreq and the number of excluded species, NR, 
as a function of time in the chemical evolution. A ‘best-fit’ solu- 
tion would be one in which D is a minimum and N, the number 
of species whose abundances are within an order of magnitude 
of those observed, a maximum. This occurs in the range of (1— 
2) x 10° yr with the formal ‘best fit’ time at 1.6 x 10° yr at which 
NR = 9 and 55 of the 134 species have calculated abundances 
that fall within an order of magnitude of their observed values 
(see Table 8). A better fit, in the sense of a lower Dyeg, can be 
found at earlier time but at the expense of removing a signif- 
icant number of species. If we seek to mimimise the number 
excluded, we again find that a time of 1.6 x 10? yr is the pre- 
ferred time. At times longer than 3 x 10? yr, NR increases and 
N decreases so that at 10° yr, D = 3.34, NR = 39 and N = 20, an 
extremely poor fit. This table also includes values in parenthe- 
ses when the reduced reaction ratefile is used. One sees that the 
exclusion of reactions makes the fit slightly worse although the 
overall differences are marginal at 1.6 x 10? yr. 

Table A.1 compares the model results for the full ratefile with 
the observed fractional abundances, including upper limits, of 
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149 species in TMC-1 at a time of 1.6 x 10? yr. It is clear from 
this table that the model, which has an initial C/O abundance = 
0.44, fails to reproduce the abundances of many species by sev- 
eral orders of magnitude, as indeed was the case for the RATE12 
models. These include species that are thought, from observa- 
tions of massive star-forming regions, to be either formed on icy 
grains or from the evaporation products of such surface reac- 
tions, species such as methanol, methyl formate, and dimethyl 
ether. 

In general, the C,,O abundances are well fit for n = 2, 3 but 
are severely over-produced for larger values as are the related 
species, HC,,O. This may be due to an overestimate of the RA 
rate coefficients for the CO-hydrocarbon ion reactions since 
these are estimated from three-body experimental measurements 
(Adams et al. 1989). 

One of the major changes between these results and RATE12 
is the huge decrease in the calculated abundance of propene, 
CH3CHCH», by ten orders of magnitude, contrary to the fact 
that its observed value in TMC-1 is 4.0 x 107? (Marcelino et al. 
2007). Previously, it was formed by the DR of C4H;* formed in 
successive radiative association reactions of H5 with C3H3* and 
C3H5*. These reactions, included in RATE12, were subsequently 
determined to have barriers in their entrance channels and there- 
fore not to proceed (Lin et al. 2013). The large decrease in the 
propene abundance has significant implications for larger hydro- 
carbons since it is a reactive species with large rate coefficients 
measured in NN reactions at low temperatures. For example, 
in terms of observed molecules in TMC-1, laboratory studies 
have shown that reactions of propene with C; form CH3;C,H 
and its isomer allenyl acetylene, HxCCCHCCH, those with 
CoH form propyne, CH3CCH, those with C4H form CH3C&H, 
those with CN form vinyl cyanide, CH,CHCN, crotononi- 
trile, CH3CHCHCN, methacrylonitrile, CH2C(CH3)CN, and 
allyl cyanide, CH2CHCH2CN (Cernicharo et al. 2022c), and 
those with CH form cyclopentadiene, c-C5He, and butadiene, 
CH5CHCHCH.». These newly formed species can then react 
to form other organic species. For example, c-C5Hg can react 
with CN to form cyanocyclopentadiene, C5H5CN, with C,H to 
form both fulvenallene, CSH4CCH», and ethynyl cyclopentadi- 
ene, CsH5CCH, and with C to form phenyl, c-CgHs5 (Cernicharo 
et al. 2022b). Similarly, allenyl acetylene, H;CCCHCCH, can 
react with CN to form cyanoacetyleneallene, H;?CCCHC3N, an 
isomer of CH3C5N (Shingledecker et al. 2021). All of these 
species are calculated to have abundances at least three orders of 
magnitude less than their observed values in TMC-1. Very few of 
these large hydrocarbons are included in the RATE12 ratefile but, 
as an example of the different predictions of the RATE22 models, 
we compare in Fig. 4 the abundances of four molecules in com- 
mon, including those of propene and butadiene, CH; CHCHCH». 
One sees clearly that this catastrophic decrease in the calculated 
propene abundance has a severe effect on those of the larger 
hydrocarbons. 


4.1.2. C-rich TMC-1 model 


One possible mechanism to enhance the formation rate of com- 
plex hydrocarbons is to adopt a model with an elemental C/O 
ratio greater than 1, simulating an environment, say, where O 
is removed from the gas in the form of H2O ice. We note 
that a selective depletion process for oxygen can also occur on 
warm dust in PDRs due to the high binding energy of O atoms 
(Esplugues et al. 2019; Róllig & Ossenkopf-Okada 2022). Here 
we adopt the values for the fractional abundance of C and O rela- 
tive to H5 of 2.8 x 107^ and 2.0 x 107^, respectively, values taken 
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Fig. 4. Time evolution of propene and species formed from it in the 
dark cloud, O-rich model for both RATE22 (solid lines) and RATE12 
(dashed lines). 


100 


D Factors 
So 


—! 


" L 
100000 1x108 
Time (yr) 


10000 


Fig. 5. Best-fit, D, and reduced best-fit, Dred, average difference 
between observed and calculated species for the full ratefile and C/O 
= 1.4. The curve labelled ‘NR’ indicates the number of species removed 
in calculating the global average, Dred, and N is the number of species 
whose abundances lie within an order of magnitude of those observed. 


Table 9. Fitting information on the comparison between observed and 
calculated abundances for the C-rich model with C/O = 1.4. 


Time (yr) D Dea NR N 

6.3x10* 1.70(1.75) 1.24(1.28) 10(10) 58(60) 
16x10? 1.68(1.67) 1.31(1.23) 910) 50652) 
5.0x10? 1.33(1.22) 1.06(0.90) 89)  67(80) 
1.0x106 1.28(1.28) 0.90(0.90) 1010) 79(80) 


Notes. Values in parentheses are for the case in which C/O = 1.1. See 
text for further details. 


from Agúndez & Wakelam (2013). Figure 5 shows the results 
and Table 9 provides global information on the fit at various 
times. In general the D and D,eq values are reduced for the C-rich 
model compared to the O-rich model. For these conditions, the 
abundances of many of the complex hydrocarbons discussed 
above now agree with observation to within an order of mag- 
nitude at 106 yr. In addition to propene, only CH3C4H, c-CsHe, 
CH35CHCHCN, and C5H35CN show differences with observation 
of more than an order of magnitude. 
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Table 10. Initial fractional abundances, f, relative to Hz for the CSE model. 


Species f Species f Species f 

He 170(-001)) CS 1.06(-05) SiO 5.02(-06) 
SiS 5.98(-05) CO 8.00(-04) C2H2 4.38(-05) 
HCN 4.00(-05) N2 4.00(-05) CH4 3.50(-06) 
NH; 6.00(-08) HO 2.55(-06) SiC, 1.87(-05) 
HCl 3.25(-07) HCP 2.50(-08) C2H4 6.85(—08) 
HF 1.70(-08) H2S 4.00(-09) SiH, 2.20(-07) 
PH; 1.00(-08) Na 1.00(-09) Mg 1.00(-05) 
Al 1.00(-07) Ar 1.00(-08) Ca 1.00(-09) 
Ti 1.00(-09) Fe 6.00(-09) 

Notes. a(b) = a x 10°. 
In this case the number of species whose abundances agree 0.0001 


to within an order of magnitude of those observed ranges from 
50 at 1.6x 10° yr to a maximum of 79 at 10° yr. In the latter 
case, only 10 species are different by more than four orders of 
magnitude and when these are removed, the average difference 
between observed and calculated values is less than a factor of 
10. If one considers the additional 15 upper limits in Table A.2, a 
further 7 species agree, that is the model gives agreement for 86 
out of a total of 149 species, or 58%. Many species showing the 
largest differences are the same in both models but in the C-rich 
case, there is a tendency to over-produce the large hydrocarbon 
chains, including C5H2, C6H2, C,O (n = 4-7), HC,O (n = 4-7) 
C, S (n = 4,5), and the cyanopolyynes with n > 7. 

Other values of C/O are used in the literature and to see the 
effect of these, we have also considered a model in which C/O = 
1.1, with an O abundance of 3.1 x 107^ (Barnum et al. 2022). The 
models show an improvement with respect to the C/O = 1.4 mod- 
els (see Table 9) and in particular at a time of 5.0 x 10? yr, when 
Dreq is about 15% smaller and N about 19% larger than the values 
for C/O = 1.4. In part, the improved fit is a result of a reduction 
in the over-production of large hydrocarbon chains discussed in 
the previous paragraph. As was also found for the O-rich case, 
the removal of the 53 endothermic reactions makes only a slight 
change in the results presented here for the full ratefile. 


4.2. IRC+10216 


We have also investigated the application of the ratefile to a 
model of the circumstellar envelope surrounding IRC+10216 in 
which the external interstellar UV field drives the chemistry. We 
assume standard physical parameters for the outflow. Motivated 
by the revised, and larger, distance to IRC+10216 suggested by 
Andriantsaralaza et al. (2022), we adopt a uniform mass-loss 
rate, M, of 3.0 x 10 Mo yr^!. We choose an expansion veloc- 
ity, ve, of 14.5 km s™!, and a power-law temperature distribution, 
T(r) = T.(r/R.) 7, with a stellar temperature T, = 2330 K and 
stellar radius R, = 5.0 x 10? cm. With a constant mass-loss rate 
and expansion velocity, the number density of H2, which is pro- 
portional to M/v,, follows an 7^? power law while the visual 
extinction to interstellar UV radiation follows an r^! distribu- 
tion. We have adopted input parent abundances from Agündez 
et al. (2020) and Van de Sande et al. (2021) augmented by metal 
atom abundances (Table 10). We calculate the chemistry from 
an initial radius of 10!^ cm out to 10/5 cm. At 10!4 cm, n(H5) = 
2.33 x10? cm?, T = 1424 K, and Ay = 1160 mag. 

Figures 6 and 7 show the fractional abundances of the 
cyanopolyynes and the Mg-terminated cyanopolyynes as a 
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Fig. 6. Fractional abundances of the cyanopolyynes for a mass-loss rate 
of 3.0 x 10? Mo yr^! as a function of radial distance from the central 
star. 
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Fig. 7. Fractional abundances of the Mg-terminated cyanopolyynes for 
a mass-loss rate of 3.0 x 10? Mọ yr! as a function of radial distance 
from the central star. 


function of radial distance. We compare calculated column 
densities with those of 61 observed daughter species, that is, 
species which are not parents, in IRC+10216 for a mass-loss 
rate of 3.0 x 107? Mo yr! in Table 11. We note, however, that 
observationally-derived column densities of the same molecule 
can differ by more than an order of magnitude as they are 
often dependent on a variety of factors: the particular excita- 
tion properties associated with observed transitions; the angular 
resolution of the telescope; whether or not the column density 
calculation includes an estimation of the size of the emit- 
ting region, often taken to be the photodissociation radius; and 
whether or not a beam-averaged column density is given. Thus 
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Table 11. Observed and calculated column densities (cm?) in IRC+10216 for a mass-loss rate of 3.0 x 10? M; yr7!. 


Species Observed Calculated Ref. Species Observed Calculated Ref. 
C 1.10(16) 1.50(16) A C5 7.90(14) 4.95(15) B 
C;H 4.60(14) 3.36115) ZE CN- 5.00(12) 2.67(10) C 
CN 1.10(15) 5.86(15) B HNC 1.10(14) 2.7714) ZE 
CH;NH 9.00(12) 4.10(10) T HCO* 3.50(12) 5.49(11) A,V 
H;CO 5.00(12) 5.71(11) U C5 1.00(15) 2.64(14) A 
C3H 5.60(15) 2.0004) ZE CH? 2.00(13) 5.54(13) A 
HCCC 3.00(12) 1.68(13) W HCCN 1.40(12) 1.89113) ZD 
SiC 6.00(13) 2.3114) N CH3CCH 1.80(13) 8.85(11) E 
CH»CN 8.40(12) 6.54(09) A,E CH3CN 4.30(13) 8.66112) ZE 
SIN 3.80(13) 3.25(12) Q CP 1.00(14) 2.27(12) L 
HCP 1.00(15) 5.86(15) L PN 1.00(13) 2.03(12) L 
H5CS 1.00(13) 1.01(12) E MgC5 1.00(12) 151(12 ZB 
C4H 5.50(15) 6.17(14) A MgC,H 2.00(12) 2.69(12) Y 
C4H7 7.00(11) 3.31(08) F C3N 3.10(14) 8.55(14) H 
C3N7 1.60(12) 1.87(11) H MgCN 74011) 6.79110) ZA 
MgNC 1.30(13) 143(1 ZA | HMgNC 6.00(11) 154(12 ZA 
C4H5 1.60(14) 6.38(14) X HC3N 1.50(15) 4.71(14) A 
C30 1.00(12) 5.80(10) S CH5CHCN 5.00(12) 1.58(12) E 
SiNC 2.00(12) 3.81(10) R CCP 1.20(12) 9.28(09) M 
CS 1.50(14) 7.81(13) I C5 1.00(14) 2.30(12) A 
CsH 2.60(14) 8.55(12) A HC4N 1.50(12) 153(13 ZD 
SiC3 4.30(12) 1.06(14) O C38 8.50(13) 7.96113) LJ 
MgC4H 2.20(13) 1.37(09) Y CH- 4.00(12) 4.22(11) E 
CoH 7.00(13) 1.03(13) E CsN 4.50(12) 5.61(13) A,D 
C;N7 3.00(12) 2.51 (12) G MgC3N 9.30(11) 3.07(10) Y 
HC3N 2.50(14) 3.76(14) A SiC4 7.00(12) 2.17(13) P 
C;H 1.50(12) 105(12 A,D CsS 3.60(13) 1.28(14) J 
CH- 2.00(12) 1.90(11) K CH 5.00(12) 3.99(12) A 
MgC6H 2.00(13) 1.44(12) Z MgC;N 4.70(12) 3.76(13) Z 
C;N^ 2.40(12) 1.6710) ZC HCN 1.00(14) 2.54(13) A 
HCN 3.00(13) 1.78(13) A 


Notes. The ‘Ref.’ column gives references for the observed values. a(b) = a x 10^. 


References. A: See references in Table 5 of Millar et al. (2000): B: Bakker et al. (1997); C: Agündez et al. (2010); D: Cernicharo et al. (2000); E: 
Agúndez et al. (2008) F: Cernicharo et al. (2007); G: Cernicharo et al. (2008); H: Thaddeus et al. (2008); I: Cernicharo et al. (1987); J: Bell et al. 
(1993); K: Remijan et al. (2007); L: Milam et al. (2008); M: Halfen et al. (2008); N: Cernicharo et al. (1989); O: Apponi et al. (1999); P: Ohishi 
et al. (1989); Q: Turner (1992); R: Guélin et al. (2004); S: Tenenbaum et al. (2006); T: Tenenbaum et al. (2010); U: Ford et al. (2004); V: Pulliam 
et al. (2011); W: Cernicharo et al. (1991a); X: Cernicharo et al. (1991b); Y: Cernicharo et al. (2019); Z: Pardo et al. (2021); ZA: Cabezas et al. 
(2013); ZB: Changala et al. (2022); ZC: Cernicharo et al. (2023); ZD: Cernicharo et al. (2004); ZE: Tuo et al. (in prep.). 


our comparison is at best only indicative of how the model 
behaves. 

We find that only one species, MgC4H, has a column den- 
sity more than 4 orders of magnitude different (less) than 
that observed. For our power-law temperature distribution, we 
find D = 1.00(0.99), Drea = 0.95(0.94), and that 37(38) of the 
61 species have calculated and observed column densities that 
differ by less than an order of magnitude, where the values 
in parentheses are from the reduced ratefile. One sees that the 
reduction in the number of reactions has a negligible effect 
on our overall fit, with negligible differences also on column 
densities (though not shown here). 

Eight species have column densities that fall more than a fac- 
tor of 100 less than those observed: the anions CN ^, C4H7, and 
C;N ;; the Mg-species MgC4H and MgC3N, CH2NH, CH2CN, 
and CCP. No species has a column density more than 100 times 
larger than observed. 

Figure 8 shows the radial distribution of the magne- 
sium polyacetylides for a mass-loss rate of 3.0 x 10? Mo,yr !. 


The abundance of the longest chain MgCgH is larger than the 
others primarily because the rate coefficient of the RA between 
Mg* and CgH; is larger than those of Mg* with the smaller 
cumulenes (see Sect. 3.24.4). 

The reason for the large discrepancy between model and 
observed abundances for MgC4H could be due to the fact that we 
do not differentiate between diacetylene, HC4H, and the cumu- 
lene H5C4. As a result, we do not include either of these as a 
parent molecule whereas in fact Fonfría et al. (2018) have found 
that diacetylene has a column density of around 1.9 x 106 cm? 
close to the star. We thus investigated a model in which we 
included diacetylene as a parent species with an initial abun- 
dance of 1.6 x 1075 but found that the column density of MgC4H 
increased only by an order of magnitude, still far too low. The 
reason appears to be that the RA rate coefficient (see Table 4) is 
simply too low. It also may be possible to synthesise the correct 
abundance if the larger protonated Mg-acetylide chains formed 
it via DR. That is, MgCeH?"/MgCsH5* + e^ — MgC4H as well 
as other products. This would imply substantial rearrangement 
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Fig. 8. Fractional abundances of the Mg-terminated polyacetylides for 
a mass-loss rate of 3.0 x 10 Mọ yr! as a function of radial distance 
from the central star. 


in the DR process and should be the subject of theoretical or 
experimental investigation. 


5. Summary 


We have outlined the major differences between our new release 
of the UMIST Database for Astrochemistry, RATE22, and the 
previous version RATEI2. The numbers of species and reactions 
have increased substantially in large part driven by the substan- 
tial number of new detections made over the past decade. Many 
Observed species are, however, excluded from the database due 
to the lack of identified gas-phase mechanisms. 

We have calculated time-dependent chemical kinetic mod- 
els for a dark cloud with physical conditions similar to those 
in TMC-1 for both O-rich and C-rich initial conditions. As 
expected, the C-rich models give a much better agreement with 
observations with 70-80 having abundances within an order of 
magnitude of the 134 observed. Our results, at least in terms 
of our ‘reduced least squares’ fit, are not affected significantly 
by the choice of C/O ratio although the lower value of 1.1 
reduces the degree of over-production of hydrocarbon chain 
molecules. We discuss reasons for the largest discrepancies and 
conclude that while grain surface formation is indicated for 
some species, the inefficient production of C3H7* and hence of 
propene, CH3CHCH,, is a key inhibitor of the production of 
larger hydrocarbons. We have also applied our network to the 
case of an outflow from a carbon-rich AGB star, finding agree- 
ment to within an order of magnitude between calculated and 
observed column densities for about 60% of the daughter species 
included in chemistry. 

In several cases, abundances are either greatly over- or under- 
produced. For dark clouds, these include species affected by 
the “propene catastrophe’ as discussed in Sect. 4.1.1. Other 
examples include the C,O and the HC,O species also dis- 
cussed in Sect. 4.1.1, and the Mg-bearing molecules observed in 
IRC+10216 (Sect. 4.2). The chemistry of all these species rely on 
two processes, radiative association and dissociative recombina- 
tion, and their calculated abundances are therefore very sensitive 
to the rate coefficients adopted for both and to branching ratios 
of the DR reactions. 

We have also considered the removal of the endothermic 
reactions identified in the KIDA database by Tinacci et al. 
(2023). We find that these make only a minor change to our 
global fits for all models considered here. 

We have not attempted to produce isotopic versions of 
RATE22, the most important of which would include 13C and 
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D. Both would need significant effort if multiply-substituted 
species were to be included. Even if one considers single iso- 
topic substitutions, there are severe uncertainties in the choice 
of additional reactions and rate coefficients. For example, during 
a formation process that involves reactants containing multiple 
carbon atoms, there is almost no information on where a °C 
atom might replace a '?C atom, nor where ?C might be retained 
in the DR of a complex hydrocarbon ion. Nearly 200 species in 
RATE22 contain 5 or more C atoms. With the exception of HC;N 
and HC7N, (Burkhardt et al. 2018), there is little observational 
data available on '?^C/"C ratios in larger molecules. Molecules 
with large numbers of carbon atoms are more likely to con- 
tain a C atom simply from a statistical viewpoint (Burkhardt 
et al. 2018), further confusing an understanding of the impact of 
5C chemistry. 

For deuterium, the situation is complicated by the need 
to include the spin chemistry of H2, its ions and deuterated 
counterparts, as well as those of multiply-deuterated species 
given that many such species are detected in interstellar clouds. 
Majumdar et al. (2017) published such a model based on a 
KIDA network that incorporated 7509 gas-phase reactions. The 
ratefile contained 111 000 reactions when deuterated. A similar 
multiplicative factor might be expected for the RATE22 ratefile. 

In addition to the gas-phase chemistry discussed here, we 
have produced a website that will include the codes needed 
to generate chemical models for both interstellar clouds and 
circumstellar envelopes. Previously, Van de Sande & Millar 
extended the RATE12 CSE models to include the effects of 
porosity, clumpy outflows and irradiation by internal UV pho- 
tons from the AGB star (Van de Sande & Millar 2019) and from 
an nearby companion star (Van de Sande & Millar 2022). We 
shall extend the capability of our CSE models by calculating 
these photorates for use with the RATE22 network. Finally, we 
note that, although there has been a decade between this and the 
previous release of the UDfA, we plan to issue updates to the 
database on an annual basis. 
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Appendix A: Comparison of observed and calculated abundances in TMC-1 


We compare observed fractional abundances in TMC-1 with calculated abundances for two models, the first for the O-rich conditions 
described in Sect. 4.1.1 (Table A.1), the second for the C-rich conditions described in Sect. 4.1.2 (Table A.2). 


Table A.1. Observed fractional abundances, relative to H5, in TMC-1(CP) and corresponding 
calculated *best fit values at 1.6 x 10? years for the O-rich model. 


Species Observed Calculated Ref. Species Observed Calculated Ref. 
OH 3.0(-07) 1.9(-08) ZG NH; 5.0(-08) 5.3(-11) ZC 
H20 <7.0(-08) 2.1(-06) ZG CH 6.5(-08) 2.5(-00 ZK 
CN 5.0(-09) 2.0(-08) ZY HCN 1.1(-08) 4.3(-07) ZW 
HNC 2.6(-08) 3.6(-07) ZW  HCNH* 1.9(-09) 2.5(-00 ZG 
H2CN 1.5(-11) 18-11) ZG HNC <3.2(-11) 9.0(-12) ZM 
CO 8.0(-05) 1.7(-04) ZW HCO 1.1(-10) 1.2(-09) U 
HCO* 9.3(-09) 5.1(-09) ZW  NoH* 4.0(-10) 4.5(-11) ZY 
CH;NH «3.6(-09) 15-11) ZU NO 3.0(-08) 2.3(-00 ZY 
H2CO 5.0(-08) 3.7(-08) F | H3COH* <3.0(-11) 31C11)) D 
CH3OH 4.8(-09) 17-15 D HS 5.0(-10) 3.5(-13) ZY 
C3H 3.0(-09) 8.5(-10) ZC C3H* 2.4(-12) 7.9(-13) G 
C3H2 1.9(-09) 1.4(-09) ZC HCCC 1.9(-10) 1.7(-11) K 
H,CCCH* 7.0(-11) 13C10 YB CH,CCH 1.0(-08) 3.2(-11) T 
HCCN 4.4(-11) 1.1(-09) H  CH3;CCH 1.2(-08) 5.6(-12) ZC 
C20 7.5(-11) 89(1) U  CH;CN 1.5(-09) 8313) O 
HCCO 7/1 (11) 7.81 | U | CH3CN 4.7(-10) 2.8(-08) O 
CH2CO 1.4(-09) 8.9(-0) D  CH3CHCH2 4.0(-09) 8.5(-22) ZZ 
HNCO 1.3(-09) 4.5(-11) D  HCNO 7.0(-12) 21-13) D 
HOCN 1.0(-11) 5.4(-13) D  CH,CO* 3.2(-11) 1.6(-11) R 
CS 3.5(-08) 6.6(-08) B  |CH3;CHO 3.5(-10) 2.1(-12) D 
CH2CHOH 2.5(-10) 3.8(-12) ZJ = H»,NCO* <4.0(-12) 61-14 D 
HCS 5.5(-10) 1.2(-11) B  HCS* 1.0(-09) 54C11)) A 
NH2CHO <5.0(-12) 0 D HCO; 4.0(-11) 5.6(-12) D 
HCOOH 1.4(-10) 6.1(-10) D  CH3OCH; 2.5(-10) 5.5(-19) ZJ 
C,H;0H 1.1(-10) 4.0(-14) ZV HCS 4.7(-09) 4.3(-11) B 
NS 1.7(-10) 177-12 E  NS* 5.2(-12) 433-12 E 
SO 1.0(-08) 3.7710 ZY C4H 8.5(-09) 2.6(-09) YA 
CjH- 2.1(-12) 1.2(-12) ZK CN 1.2(-09) 16(10 YA 
CN- 6.4(-12) 3.1(-13) YA C4H2 3.3(-10) 23(C1) K 
HC3N 2.3(-08) 6.800 C HNC; 5.2(-11) 6.4(12 C 
HCCNC 3.0(-10) 85(10 C  HC4NH* 1.0(-10) 1.2(-11) J 
HCCNCH* 3.0(-12) 19-12 W C30 1.2(-10) 8.011) U 
CH5CHCCH 1.2(-09) 9.1(-13) H  HC3O 1.3(-11) 5.8(-10) U 
HC30* 2.1(-11) 8.2(-12) D |CH;CHCN 6.5(-10) 5.1(-14) ZC 
NCCNH* 8.6(-12) 2.7(--14) Q  HCCCHO 1.5(-10) 10(-08 F 
HCOCN 3.5(-11) 2.4-10 F | C9H5CN 1.1(-11) 1.6(-10) H 
CH,CHCHO 2.2(-11) 15(-13 ZJ CS 5.5(-09) 1.6(-11) <A 
HCCS 6.8(-11) 7.7(-12) B  HC,S* 1.1(-10) 7.0(-13) M 
NCS 7.8(-11) 1.7(-12) B  CH3COCH; 1.4(-11) 16(16) ZV 
H;CCS 1.8(-10) L8(-1) B  HSCN 5.8(-11) 2.4-13 B 
HNCS 3.8(-11) 4.6(-13) B OCS <1.8(-09) 1.2(-10) ZC 
HCOOCH; 1.1(-10) 15(-160 ZJ  CsH 1.3(-10) 3.3(-10) ZE 
CsH* 8.8(-12) 18-11) G  c-CsH 9.0(-12) 2.1(-12) ZE 
C5H5 1.4(-12) 4.0(-11) K  c-C4HCCH 3.1(-11) 17(-12) ZB 
HC4N 3.7(-11) 11-1) H C40 «9.0(-12) L4(-089 U 
SO» 3.0(-10) 1.1(-13) ZX  CH3C4H 1.3(-09) L4C13 S 
H?CCCHCCH 1.2(-09) 1.4(-13) S | CH?C;N 1.6(-11) 91(-1) N 
HC4O «9.0(-12) 79(-00 U  CH3C3N 1.7(-10) 1.6(-10) P 
HCCCH;CN 2.8(-10) 4.6(-14) P | H;CCCHCN 2.7(-10) 5.3013) P 
C5Hg 1.2(-09) 11-20 ZB  CH3CHCHCN 1.8(-11) 6.2(-24) ZO 
C38 1.3(-09) 17(-00 A  HCsSS «2.9(-12) 3.2(-10) B 
HC3S* 4.0(-11) 11-11) B  H3CCCS 7.8(-12) 62(-1) B 
HCSC,H 3.2(-11) 44-11) F | HCSCN 1.3(-10) 8.9(-13) F 
CoH 4.8(-10) 1510 YA C6H7 1.6(-11) 3.8(-12) ZK 
C6H2 8.0(-12) 1.1-1) K  CsN 4.7(-11) 233-1) YA 
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Species Observed Calculated Ref. Species Observed Calculated Ref. 
CsN- 8.8(-12) 1.2(-12) YA HCsN 1.8(-08) 4.7(-09) C 

HC4NC 3.0(-11) 18C10 | C  HCSNH* 7.5(-11) 2.9(-12) J 

C6H4 5.0(-11) 2.8(-12) Y C50 1.5(-12) 1.4(-08) D 

HC5O 1.4(-10) 6.4.00 U  NC4NH* 1.1(-12) 1.6(-12) ZR 
H2CCHC3N 2.0(-11) 1.2(--14) ZI  HCCCHCHCN 3.0(-11) 2.2(-14) ZI 
C4S 3.8(-12) 8.2(-10) B  HC4S 9.5(-12) 5.7(-13) ZQ 
C;H 6.5(-12) 5.1(-11) ZK  CgO «L1(-11) 8.8-00 U 

CH3C6H 7.0(-11) 2.2(-14) ZF  H;CCCHC4H 2.2(-10) 1.4(-11) ZF 
CH3C5N 9.5(-12) 2.3(-11) ZF H,CCCHC3N 1.2(-11) 1.3(-15) ZF 
HCgO <1.8(-11) 4.4(-09) U  CsHsCCH 3.4(-10) 9.6(-14) ZP 
CsH4yCCH2 2.7(-10) 1.6(-13) ZK  CsHsCN 1.0(-10) 15(-20 Z 

CsS 5.0(-12) 1.8(-12) B CH 4.6(-11) 4.3(-11) ZH 
CH- 2.1(-12) 1.1(-12) ZH CN% 5.0(-12) 2.9(-13) ZT 
HCN 6.4(-09) 2.4(-10) C  HC;NH* 5.5(-12) 69(-13 L 

C;O <2.6(-12) 6.6(-09) ZN  HC;O 6.5(-11) 5.3(-00 U 

C;H5CCH 2.5(-10) 6.9(-13) ZP CeHsCN 1.6(-10) 7.6(-13) ZA 
CoH <3.5(-12) 12-11) ZK CHCH <9.8(-10) 5.9(-21) V 

CH3C7N 8.6(-12) 9.7(-13) V | CjgH- 4.0(-11) 2.8(-13) ZS 
CioH 2.0(-11) 14-11) ZS  HCoN 1.1(-09) 9.9(-11) ZC 
HCN 1.0(-10) 2.2(-12) ZD 


Notes. The ‘Ref.’ column gives references for the observed values. a(b) = a x 10°. 


References: A: Cernicharo et al. (2021j); B: Cernicharo et al. (2021g); C: Cernicharo et al. (2020a); D: Cernicharo et al. (2020b); E: Cernicharo 
et al. (2018); F: Cernicharo et al. (20211); G: Cernicharo et al. (2022a); H: Cernicharo et al. (2021a); I: Cernicharo et al. (2020c); J: Marcelino 
et al. (2020); K: Cabezas et al. (2021d); L: Cabezas et al. (2022c); M: Cabezas et al. (2022b); N: Cabezas et al. (2021b); O: Cabezas et al. (2021c); 
P: Marcelino et al. (2021): Q: Agündez et al. (2015a); R: Cernicharo et al. (2021h); S: Cernicharo et al. (2021f); T: Agündez et al. (2022b); U: 
Cernicharo et al. (2021b); V: Siebert et al. (2022a); W: Agündez et al. (2022a); Y: Cernicharo et al. (2021d); Z: Lee et al. (2021a); ZA: Burkhardt 
et al. (2021b); ZB: Cernicharo et al. (2021c); ZC: Gratier et al. (2016); ZD: Loomis et al. (2021); ZE: Cabezas et al. (2022a); ZF: Fuentetaja et al. 
(2022b); ZG: Agúndez & Wakelam (2013); ZH: Brünken et al. (2007); ZI: Lee et al. (2021b); ZJ: Agúndez et al. (2021b); ZK: Cernicharo et al. 
(2022b); ZL: Loison et al. (2016); ZM: Cabezas et al. (2021a); ZN: Cordiner et al. (2017); ZO: Cernicharo et al. (2022c); ZP: Cernicharo et al. 
(2021e); ZQ: Fuentetaja et al. (2022a); ZR: Agündez et al. (2023a); ZS: Remijan et al. (2023); ZT: Cernicharo et al. (2023); ZU: Kalenskii et al. 
(2004); ZV: Agündez et al. (2023b); ZW: Pratap et al. (1997); ZX: Cernicharo et al. (2011); ZY: McElroy et al. (2013); ZZ: Marcelino et al. (2007); 


YA: Agündez et al. (2023c); YB: Silva et al. (2023). 
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Table A.2. Observed fractional abundances, relative to H2, in TMC-1(CP) and corresponding 
calculated ‘best fit’ values at 1.0 x 10° years for the C-rich (C/O = 1.4) model. 


Species Observed Calculated Ref. Species Observed Calculated Ref. 
OH 3.0(-07) 3.6(-08) ZG NH; 5.0(-08) 4.7(-08) ZC 
H20 <7.0(-08) 9.4(-08) ZG CH 6.5(-08) 6.1(-09) ZK 
CN 5.0(-09) 7.2(-08) ZY HCN 1.1(-08) 1.4(-07) ZW 
HNC 2.6(-08) 18007 ZW HCNH* 1.9(-09) 8.9(-10 | ZG 
H2CN 1.5(-11) LOCI) ZG HNC <3.2(-11) 4.2(-12) ZM 
CO 1.7(-04) 1.9(-04) ZW HCO 1.1(-10) 5.2(-11) U 
HCO* 9.3(-09) 5.5(-00 ZW  NoH* 4.0(-10) 7.0(-10) ZY 
CH2NH <3.6(-09) 3.5(-00 ZU NO 3.0(-08) 5.3(-08) ZY 
H,CO 5.0(-08) 5.7(-00 F  H;COH* <3.0(-11) 4.1(-12 D 
CH30H 4.8(-09) 2.7(--13) D HS 5.0(-10) 5.6(-10) ZY 
C3H 3.0(-09) 3.3(-09) ZC C3Ht 2.4(-12) 4.0(-14) G 
C3H5 1.9(-09) 8.3(-09) ZC HCCC 1.9(-10) 1.4(-11) K 
H2CCCH* 7.0(-11) 1.0(-10) YB  CH2CCH 1.0(-08) 3.0(-11) T 
HCCN 4.4(-11) 7.5(-10) H  CH3CCH 1.2(-08) 1.8(-09) ZC 
C20 7.5(-11) 1.1(-09 U  CH5;CN 1.5(-09) 11(-1) O 
HCCO 7.1(-11) 7.5(-13) U | CH3CN 4.7(-10) 4.6(-10) O 
CH2CO 1.4(-09) 7.2(-10) D  CH;CHCH, 4.0(-09) 3.4(-16) ZZ 
HNCO 1.3(-09) 33-1) D | HCNO 7.0(-12) 5.3(-13) D 
HOCN 1.0(-11) 5.60212 D  CH3CO* 32(-11) 9.201135 R 
CS 3.5(-08) 110-07) B  CH3CHO 3.5(-10) 5.1012 D 
CH,CHOH 2.5(-10) 1.4(-11) ZJ H3NCO* <4.0(-12) 3.5(-12 D 
HCS 5.5(-10) 2.3(-12) B  HCS* 1.0(-09) 6.5(-11) A 
NH2CHO <5.0(-12) 0 D HCO; 4.0(-11) 4.3(-12) D 
HCOOH 1.4(-10) 2.5(-11) D CHOCH; 2.5(-10) 9.2(-19) ZJ 
C2H5OH 1.1(-10) 7.0(-13) ZV  H3CS 4.7(-09) 12(08 B 
NS 1.7(-10) 5.0-10 | E  NS* 5.2(-12) 73(12) E 
SO 1.0(-08) 3.9(-08) ZY C,H 8.5(-09) 9.4(-09) YA 
C4H- 2.1(-12) 6.3(-13) ZK CN 1.2(-09) 2.1(-09) YA 
CN- 6.4(-12) 6.7(-13) YA C4H2 3.3(-10) 2.3(-09) K 
HC3N 2.3(-08) 1.1(-08) C HNC; 5.2(-11) 15-11) C 
HCCNC 3.0(-10) 82C10 C  HC3NH* 1.0(-10) 9.8(-12) J 
HCCNCH* 3.0(-12) L6-12 W C30 1.2(-10) 8.9(-0 U 
CH,CHCCH 1.2(-09) 6.6(-10) H  HC3O 1.3(-11) 53(11)) U 
HC30* 2.1(-11) 2.8012 D | CH;?CHCN 6.5(-10) 7.3(-10) ZC 
NCCNH* 8.6(-12) 14-11) Q  HCCCHO 1.5(-10) 2.5(-10) F 
HCOCN 3.5(-11) 6.8(-11) F = C,HsCN 1.1(-11) 3.0(-11) H 
CH5CHCHO 2.2(-11) 16(-13 ZJ GS 5.5(-09) 51000 A 
HCCS 6.8(-11) 11010 B  HC,S* 1.1(-10) 1.3(-11) M 
NCS 7.8(-11) 7.2(-11) B | CH3COCH; 1.4(-11) 2.8(-17) ZV 
H2CCS 1.8(-10) 49-11) B  HSCN 5.8(-11) 19-1) B 
HNCS 3.8(-11) 33C-1) B OCS <1.8(-09) 2.7(-08) ZC 
HCOOCH; 1.1(-10) 8.0(-19 ZJ CsH 1.3(-10) 5.3(-00) ZE 
CsH* 8.8(-12) 18-11) G  c-Cs5sH 9.0(-12) 5.3(-00) ZE 
C5H5 1.4(-12) 1.4-0) K  c-C;HCCH 3.1(-11) 4.8(-11) ZB 
HC4N 3.7(-11) 1.3(-09) H C40 «9.0(-12) 22(-00 U 
SO» 3.0(-10) 2.8-00 ZX CH3C4H 1.3(-09) 61(-11) S 
H2CCCHCCH 1.2(-09) 14-10) S | CH2C;N 1.6(-11) 18-11) N 
HC4O «9.0(-12) 3.4-00 U  CH3C;N 1.7(-10) 12(-00 P 
HCCCH;CN 2.8(-10) 8.2(-11) P | H5CCCHCN 2.7(-10) 8.8(-10 P 
C5H6 1.2(-09) 6.8(-15) ZB  CH3CHCHCN 1.8(-11) 2.4(-17) ZO 
C38 1.3(-09) 9.4-00 A  HC3SS «2 A11) 14-11) B 
HC3S* 4.0(-11) 47-11) B  H3;CCCS 7.8(-12) 47(M) B 
HCSC,H 3.2(-11) 2.3(-10) F = HCSCN 1.3(-10) 7.0(-10) F 
CoH 4.8(-10) 6.3(-09) YA CH7 1.6(-11) 2.6(-10) ZK 
C6H2 8.0(-12) 5.007 K  CsSN 4.7(-11) 9.8(-11) YA 
CsN~ 8.8(-12) 7.2(-12) YA | HC3N 1.8(-08) 2.5000 C 
HC4NC 3.0(-11) 4.8(-10) C  HCjNH* 7.5(-11) 1.0(-11) J 
C6H4 5.0(-11) 7.3(-10) Y C50 1.5(-12) 1.3(-09) D 
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Species Observed Calculated Ref. Species Observed Calculated Ref. 
HC;0 1.4(-10) 17709 U  NC4NH* 1.1(-12) 5.2(-11) ZR 
H2,CCHC3N 2.0(-11) 7.2(-11) ZI HCCCHCHCN 3.0(-11) 1.4(-10) ZI 
C48 3.8(-12) 5.3-00 B  HC,S 9.5(-12) 11(-00 ZQ 
C;H 6.5(-12) 2.5(-00 ZK C,.O «L1(-11) 1.4(-08) U 
CH3C6H 7.0(-11) 8.0(-11) ZF  H3CCCHC4H 2.2(-10) 1.1(-10) ZF 
CH3C5N 9.5(-12) 1.1(-11) ZF H,CCCHC3N 1.2(-11) 8.9(-12) ZF 
HC,O <1.8(-11) 6.2(-09) U  CsHsCCH 3.4(-10) 4.1(-11) ZP 
CsH4CCH2 2.7(-10) 2.7(-10) ZK  CsHsCN 1.0(-10) 2.5(-15) Z 
C5S 5.0(-12) 2.4.00 B CH 4.6(-11) 7.2(-09) ZH 
CH- 2.1(-12) 1.7(-10) ZH CN 5.0(-12) 2.0010) ZT 
HCN 6.4(-09) 4.007 C  HC;NH* 5.5(-12) 6.5(-10) L 
C70 <2.6(-12) 1.3(-07) ZN  HC;O 6.5(-11) 12(07 U 
C6H5CCH 2.5(-10) 3.4(-10) ZP  CęH5CN 1.6(-10) 2.0010) ZA 
CoH <3.5(-12) 3.4(-09) ZK  CH3C;H <9.8(-10) 3.0012 V 
CH3C7N 8.6(- 12) 117-10) V  CygH- 4.0(-11) 1.1(-10) ZS 
CioH 2.0(-11) 4.5(-09) ZS | HCo9N 1.1(-09) 5.0(-07 ZC 
HCjN 1.0(-10) 17(07 ZD 


Notes. The ‘Ref.’ column gives references for the observed values. a(b) = a x 10^. 


References: A: Cernicharo et al. (2021j); B: Cernicharo et al. (2021g); C: Cernicharo et al. (20202); D: Cernicharo et al. (2020b); E: Cernicharo 
et al. (2018); F: Cernicharo et al. (20211); G: Cernicharo et al. (2022a); H: Cernicharo et al. (2021a); I: Cernicharo et al. (2020c); J: Marcelino 
et al. (2020); K: Cabezas et al. (2021d); L: Cabezas et al. (2022c); M: Cabezas et al. (2022b); N: Cabezas et al. (2021b); O: Cabezas et al. (2021c); 
P: Marcelino et al. (2021): Q: Agúndez et al. (20152); R: Cernicharo et al. (2021h); S: Cernicharo et al. (2021f); T: Agúndez et al. (2022b); U: 
Cernicharo et al. (2021b); V: Siebert et al. (2022a); W: Agúndez et al. (2022a); Y: Cernicharo et al. (2021d); Z: Lee et al. (2021a); ZA: Burkhardt 
et al. (2021b); ZB: Cernicharo et al. (2021c); ZC: Gratier et al. (2016); ZD: Loomis et al. (2021); ZE: Cabezas et al. (2022a); ZF: Fuentetaja et al. 
(2022b); ZG: Agúndez & Wakelam (2013); ZH: Brünken et al. (2007); ZI: Lee et al. (2021b); ZJ: Agúndez et al. (2021b); ZK: Cernicharo et al. 
(2022b); ZL: Loison et al. (2016); ZM: Cabezas et al. (2021a); ZN: Cordiner et al. (2017); ZO: Cernicharo et al. (2022c); ZP: Cernicharo et al. 
(2021e); ZQ: Fuentetaja et al. (2022a); ZR: Agúndez et al. (20232); ZS: Remijan et al. (2023); ZT: Cernicharo et al. (2023); ZU: Kalenskii et al. 
(2004); ZV: Agúndez et al. (2023b); ZW: Pratap et al. (1997); ZX: Cernicharo et al. (2011); ZY: McElroy et al. (2013); ZZ: Marcelino et al. (2007); 


YA: Agündez et al. (2023c); YB: Silva et al. (2023). 
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